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1. Introduction 

The Mumford-Shah functional was introduced in [36] and [37] within the con- 
text of a variational approach to image segmentation problems (cf. [37] and [31], 
Chapter 4). In dimension n it can be written as follows 

F(u):= I \Vu\ 2 dx + aH^iSu) + f3 f(u-g) 2 dx , (1.1) 

where Cl is a bounded regular domain in R n , g : Q — > [0, 1] is a given function 
(input grey level), a and (3 are positive (tuning) parameters, 7^ n_1 is the (n — 1)- 
dimensional Hausdorff measure (that is, the usual (n— 1) -dimensional area in case 
of subsets of regular hypersurfaces, the length in the most relevant case n = 2). 
The unknown function u : £1 — > R is regular (say, of class C 1 ) out of a closed 
singular set Su, whose shape and location are not prescribed. Thus minimizing F 
means optimizing the function u and the singular set Su. Indeed, in the original 
formulation only the planar case n = 2 is considered, and the singular set is 
explicitly viewed as an independent variable, u being smooth on the complement 
of this set. 

While existence results for minimizers of F in dimension two could be proved 
within the original framework (cf. [37], [18], or [31], Chapter 15), in arbitrary 
dimension they were first obtained by a different approach outlined by E. De 
Giorgi (cf. [20]). More precisely, F can be defined for every function u in the 
space SBV(Q) of special functions with bounded variation (see Section 2 for more 
details, or [6], Chapter 4; this space can be regarded as a sort of the completion of 
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piecewise regular functions), upon which it is lower semicontinuous and coercive 
with respect to the L 1 topology (cf. [3], or [6, Chapter 6). This immediately 
yields minimizers of F within this space, which can be proved to be minimizers 
of the original functional by suitable regularity theorems (cf. [21], see also [6], 
Chapter 6, for further regularity results). Furthermore, the lower semicontinuity 
and coercivity results in [3] apply to a larger class of functionals coupling bulk 
and surface contributions, among which F should be regarded as a prototypical 
example. Consequently, the SBV setting has been used to model a wide range of 
problems, from image segmentation, to fractures in brittle materials, to nematic 
liquid crystals (see [6], Section 4.6, for a survey). 

On a mathematical level, one of the most relevant features of the Mumford- 
Shah functional is a deep lack of convexity. Hence, not only minimizers may be 
not unique, but "identifying" them is by no means an easy task, also in terms of 
efficient algorithms. Clearly, every minimizer u of F must satisfy certain equilib- 
rium conditions - an Euler-Lagrange equation of a sort - which can be obtained 
by considering different types of infinitesimal variations (see for instance [37] or 
[6], Section 7.4). Among these we mention the following: u satisfies Au = j3(u — g) 
in the complement of the singular set Su, the normal derivative of u on Su must 
vanish (where Su is a regular surface), while the mean curvature of Su multiplied 
by a is equal to the difference of the energy densities \ Vu\ 2 + (3{u — g) 2 on the two 
sides of Su (the first two correspond respectively to the Euler-Lagrange equation 
for the convex functional f^ K [\ Vu\ 2 + (3(u — g) 2 ] dx and the associated Neumann 
boundary condition on dK; together they are equivalent to minimality with pre- 
scribed singular set Su = K). However, due to the lack of convexity of F, these 
conditions do not imply minimality - not even local minimality. 

In this paper we propose a sufficient condition for minimality, and describe 
some applications. 

We remark that the problem of finding sufficient conditions for minimality 
already makes sense for a simplified version of F which occurs in the theory of 
interior regularity for minimizers of F, and is obtained by dropping the lower 
order term in (1.1) and setting for simplicity a = 1, that is, 



For the time being we focus on minimizers of Fq with prescribed boundary values 
(in short, Dirichlet minimizers), and describe the basic idea of this paper, without 
dwelling on details. More precise definitions and accurate statements will be given 
in Section 3. 

Assuming that u and Su are sufficiently regular, let u + and u~ denote the 
limits of u on the two sides of Su, so that u + > u~ , and let u u be the unit normal 
to Su pointing from the side of u~ to that of u + ; the complete graph of u is the 
boundary of the subgraph of u (the set of all points (x, t) € Q x R such that 
t < u(x)), oriented by the inner normal v Vu . Thus Tu consists of the union of 
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the usual graph of u and an additional part given by all segments with endpoints 
(x,u~(x)) and (x,u + (x)), with x ranging in Su; see Figure 1 (cf. also Remark 
2.11). 
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Figure 1 



We consider now the vectorfields cf) = (</> x ,(/>*) on x R such that, for every 
function u, F (u) is larger than or equal to the flux of <p through Fu, that is 

F (u) > [ <t>-v Vu dH n . (1.3) 

JYu 

Taking into account (1.2) and the fact that the flux of 4> through Fu is given by 
(cf. Figure 1) 

/ [(f> x (x,u) ■ V«-^(i,m)] dx+ [ \f 4> x {x,t)dt\-u u dH n - 1 , (1.4) 
Jn J Su *-J u - J 

we see that inequality (1.3) is satisfied if 

4> x (x, u) • Vu - ^(XjU) < |Vu| 2 for every x G il\Su 

and 



J u 



: (x,t) dt 



< 1 for every x G Su. 



The first inequality holds for every u if <j> satisfies (f) x (x,t) ■ £ — </>*(x,t) < |£| 2 for 
every point (x, t) and for every vector £, which can be equivalently restated as 

(a) \(j) x (x,t)\ 2 < 44*(x,t) for x G fi, t G R, 

while the second inequality is satisfied if 

ft 2 



(b) />(x,t) 



dt 



< 1 for x G fi, ti,t 2 G R. 



Moreover, one easily checks that equality holds in (1.3) for a particular u if (and 
only if) 

(a') (j) x (x,u(x)) = 2Vu(x) and <f>\x,u(x)) = |Vu(x)| 2 for x G Q \ Su, 



pU~t~ (x) 

(b') / (f) x (x,t) dt = v u (x) for x G Su. 
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Let now be given a function u, and assume that there exists a vectorfield 4> 
which is divergence-free and satisfies assumptions (a), (b), (a'), and (b') above. 
Then for every function v which agrees with u on the boundary of Q we have 



where the first equality follows from the divergence theorem, since (ft is divergence- 
free and Tu and Tv have the same boundary. Hence the existence of such a 
vectorfield <f> implies the minimality of u. In Section 3 we give a more precise 
version of this result, and extend it to minimizers of F and other functionals 
(cf. Theorems 3.4 and 3.10). In particular, for minimal partitions we recover 
the principle of paired calibrations introduced in [34] and [11] (in fact, a slight 
generalization of it - cf. Theorem 3.11 and Remark 3.12). 

We call 4> a calibration (for u) by analogy with the corresponding theory for 
minimal surfaces: in that setting, a vectorfield 4> is said to calibrate an oriented 
hypersurface S (with boundary) if it agrees on S with the normal vectorfield, is 
divergence- free, and satisfies \<j)\ < 1 everywhere; the existence of a calibration 
implies that S minimizes the area among all oriented hypersurfaces with the 
same boundary, and the proof is just one line, as (1.5) above (we refer the reader 
to [32] for detailed references and a review of many results). Calibrations have 
also been defined for general integrals on oriented surfaces (or currents) of any 
dimension and codimension (cf. [24], notice that in this general framework they 
are closed differential form rather than divergence-free vectorfields). Moreover, 
since every variational integral of the form f(x, u, Vu) dx with u : Q — > R fe and 
/ polyconvex in Vu can be re-written as a one-homogeneous convex integral on 
the graph of u, the theory of calibrations can be adapted also to these problems, 
yielding the equivalent of null-lagrangians. 

However, this geometric interpretation does not cover free-discontinuity prob- 
lems, because functionals of Mumford-Shah type cannot be written as integral 
over the (complete) graph of u. Indeed, the novelty of our approach consists in 
introducing suitable non-local constraints (namely, condition (b)) to define the 
class of admissible calibrations. 

Once the general principle of calibrations is given, two fundamental questions 
arise: does every minimizer admit a calibration? And how can we recover it? In 
other words, we have given a sufficient condition for minimality, but we do not 
know if it is actually fulfilled by any minimizer, and, which is even more relevant 
in applications, we do not know how to verify it, that is, how to construct a 
calibration. 

To discuss the first question, we begin by recalling the situation of minimal 
surfaces. The basic idea behind calibrations is that the area functional admits a 
natural extension from regular oriented surfaces (with fixed boundary) to all nor- 
mal currents (with the same boundary), which is nothing else but the mass. But 
now we have a convex functional on an affine space, and hence minimum points 
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are exactly characterized by the fact that the subdifferential of the functional 
contains the zero element; in the specific case, this turns out to be equivalent to 
the existence of a calibration with Borel coefficients (cf. [24], Proposition 4.10(3)). 
Thus a given surface admits a calibration of a sort if and only if it minimizes the 
area (the mass) among all normal currents with same boundary. However, in 
codimension (and dimension) larger than one, minimizing the area in the class of 
surfaces, or even integral currents, does not necessarily imply minimizing the area 
in the class of normal currents. In other words, the infima of the area functional 
on the two classes may not coincide: counterexamples were given, for instance, in 
[33] and [40]. 

The situation of the Mumford-Shah functional is somehow similar. Let Go be 
the class of all vectorfields cj> on Q. x R which satisfy assumptions (a) and (b) above; 
since for every function u we can find <fi 6 Go which also satisfies assumptions (a') 
and (b'), we have the identity 

F Q (u) = sup f • v Tu dH n = sup f <P- Dl u , (1.6) 

where Dl u is the (distributional) derivative of the the characteristic function 1„ of 
the subgraph of u. Thus F (u) = G(l u ), where G(v) is defined as the supremum 
of Jq xR </> • Dv over all (j) £ Go f° r every function v G BV\ oc {Vi x R). Since G 
is a convex functional on the affine space X u of all functions v which agree with 
l u on the boundary, l u minimizes G if and only if the subdifferential of G at l u 
contains the zero element, condition which is roughly equivalent to the existence 
of a calibration for u. It was proved in [15] for the one-dimensional case n = 1 
that the function l u minimizes G whenever u minimizes Fq; in other words, the 
infimum of G on X u agrees with the infimum of F (among all functions with same 
boundary values as u) , thus proving that every minimizer of Fq can be calibrated 
(in some sense). It is not known if the same result holds in higher dimension. 

The discussion of this point may benefit from being set in a more abstract 
framework. We can indeed summarize the calibration method as follows: given 
a nonconvex function(al) F on a certain space X, we embed X into a convex 
space X so that F agrees with a convex function F on X. Thus the minimality 
on a point x G X with respect to F can be characterized via the subdifferential 
of F at x; this yields a sufficient condition for the minimality with respect to F, 
conditions that is also necessary when (and only when) the infima of F on X and 
X agree. As pointed out in [10], it is always possible to construct the convex 
space X and the function F so that the two infima agree, but the problem is that 
checking the condition of minimality derived from this abstract construction may 
be neither simpler nor more feasible than a direct verification of minimality. 

Passing to the second question mentioned above, we notice that an abstract 
result which guarantees the existence of calibrations would nevertheless provide 
no solution to the problem of construction. As a matter of fact, we do not know 
of any general method of construction, not even for minimal surfaces. Instead, we 
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have collected in Sections 4 and 5 many examples of calibrations for F and F, 
and gathered some helpful remarks and observations. Despite a lack of a general 
recipe, the calibration method provides short and easy proofs of some natural 
minimality result; among these we recall the following ones: 

(1) every harmonic function is a Dirichlet minimizer of |Vu| 2 dx, and is 
also a Dirichlet minimizer of F when the gradient is sufficiently small 
(Paragraph 4.10); 

(2) a function which is constant on each element of a minimal partition of 
the domain is a Dirichlet minimizer of Fq when the values are sufficiently 
far apart from each other; in particular this applies to the so-called triple 
junction (Paragraphs 4.14 and 4.16); 

(3) every solution of the equation — Au + (3(u — g) = with Neumann bound- 
ary conditions is a minimizer of Vii| 2 + (3(u — g) 2 ] dx, and also of F 
for large f3 (Paragraph 5.4); 

(4) if g is the characteristic function of a regular set, then u := g minimizes 
F for large (3 (Paragraph 5.5). 

Notice that (3) and (4) give a strong indication that for initial data with smooth 
singular sets the gradient flow associated with F in the L 2 metric (which can be 
defined via time discretization, cf. [27]) leaves the singular set still, at least for 
small times, and agrees with the heat flow elsewhere (cf. [26] and Remarks 5.12, 
5.15, 5.16). 

We conclude this introduction with a few remarks. A long-standing conjecture 
on the Mumford-Shah functional is that the "cracktip" , namely the function u on 
the plane given in polar coordinates by u := \f2ppK sin(#/2), with — 7r < 6 < tt, 
(or, equivalently, the imaginary part of \J~2zpn with a cut along the negative real 
axis) is a Dirichlet minimizer of the homogeneous Mumford-Shah functional Fq 
on every bounded open subset of R 2 . This conjecture has been recently proved 
in [9], but so far no calibration has been found for this minimizer. 

As shown in [12] in the case of minimal partitions, the calibration method can 
be numerically implemented to get rigorous lower bounds for the value of minima, 
to disprove the minimality of a given configuration, or viceversa to get an idea of 
what the calibration (if it exists) should look like. It is not clear if an efficient 
numerical implementation for the Mumford-Shah functional is feasible. 

The theory presented in this paper is limited to scalar functions. To develop a 
similar theory for functionals on Revalued maps u, one should replace divergence- 
free vectorfields by closed n- forms on Q x R fc , which act on the graphs of the maps 
u, viewed as (suitably defined) n-surfaces in Q, x R fe . 

This paper is organized as follows. In Section 2 we recall the basic notation 
about finite perimeter sets and the space SBV, which is indeed the natural setting 
for our theory (however, under most regards the unfamiliar reader can just replace 
the word SBV with "smooth out of a piecewise smooth singular set", and "finite 
perimeter" with "piecewise smooth boundary" ) . In Section 3 we expand the idea 
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outlined above and develop the theory of calibrations for minimizers of F and F, 
with or without prescribed boundary values. Then we discuss the extension to 
more general functionals, and the connection with paired calibrations. Sections 4 
and 5 are devoted to examples and applications. Finally, the appendix contains 
the proofs of some technical results stated in Section 2. 

Some of the results contained in this paper were announced in [2] and [16]. 
Further applications can be found in [17] and [30]. 

Acknowledgements. - The first and third authors have been partially sup- 
ported by MURST through the projects "Equazioni Differenziali e Calcolo delle 
Variazioni (1997)" and "Calcolo delle Variazioni (2000)". This research was initi- 
ated while the first author was visiting the University of Toulon, and subsequently 
developed during a stay at the Max Planck Institute for Mathematics in the Sci- 
ences in Leipzig. Several people contributed, with discussions and remarks, to 
the final shape of this paper; among them, we would like to thank in particular 
Antonin Chambolle and Massimo Gobbino. 

2. Notation and preliminary results 

Throughout this paper, sets and functions are always assumed to be Borel 
measurable, and we do not identify functions which agree almost everywhere. A 
vectorfield on a subset E of R n is any map from E into R n . The divergence of 
a vectorfield is always intended in the sense of distributions (in the interior of its 
domain) ; in particular we say that 4> is divergence- free to mean that its divergence 
vanishes. 

The letter denotes a (possibly unbounded) open subset of R n , S n_1 is the 
unit sphere in R™, TL k stands for the /c-dimensional Hausdorff measure (which 
agrees with the usual /c-dimensional volume on every regular surface of dimension 
k), and C n is the n-dimensional Lebesgue measure; when the integration is done 
with respect to Lebesgue measure, we always write dx instead of dC n . The re- 
striction of any Borel measure \i to a set E is denoted by [i\_E, while g ■ /j, is the 
(vector) measure canonically associated with any ^-summable (vector) function 
g. Thus fiL£ = 1e ■ \i where 1e is the characteristic function of E, namely the 
function which is equal to 1 on E and to elsewhere. || • || p denotes the norm in 
the space L p . 

A (vector) function / on R n has approximate limit a at x, and we write 
aplim/(y) = a, if 



The same definition applies to functions defined on a subset E of R n , provided 
that B(x, r) is replaced by E n B(x, r). Notice that this definition slightly differs 
from the more usual one, which is expressed in term of the density at x of the 
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pre-images of neighbourhoods of a (cf. [22], Section 1.7.2, and [6], Definition 3.63); 
however, these two notions agree for locally bounded functions. 

We recall now some notation and basic facts about finite perimeter sets, BV 
and SBV functions; for a more precise definitions and a detailed account of the 
results we refer to [6], Chapters 3 and 4 (for the theory of BV functions, see also 
[22], Chapter 5). 

A real function u on an open set £1 C R n has bounded variation in f2, and we 
write u 6 BV(Cl), if it belongs to L 1 ^) and the distributional derivative Du is 
(represented by) a bounded vector measure on £1. The integral of an R n -valued 
function / with respect to the measure derivative Du is denoted by J n f- Du. 
Every u G BV(Q) is almost everywhere different iable in the approximate (or 
measure theoretic) sense, and the corresponding approximate gradient Vm agrees 
with the density of the measure Du with respect to Lebesgue measure. If O has 
Lipschitz boundary, then u admits a trace on the boundary (in the approximate 
sense), which we still denote by u, and which belongs to Ll oc (dQ,TC n ~ 1 ). 

The singular set Su is the set of all points where u admits no approximate 
limit; Su has Hausdorff dimension (less than or) equal to n — 1, and more precisely 
it is rectifiable (of dimension n— 1), which means that it can be covered, up to an 
H n_1 -negligible subset, by countably many hypersurfaces of class C 1 (these sets 
are sometimes called "countably (W n_1 , re— l)-rectifiable" , cf. [23], Section 3.2.14). 
We recall that, for every rectifiable set, the approximate unit normal, and the 
corresponding approximate tangent hyperplane, are well-defined at H n_1 -almost 
every point (and do not depend on the choice of the covering). The approximate 
unit normal to Su at x is denoted by ^„(x). At H ra_1 -almost every x G Su there 
exist the approximate limits u + (x) and u~ (x) of u on the two sides of Su (more 
precisely, the approximate limits of the restriction of u to the two half-spaces 
defined by the approximate tangent hyperplane at x), and we arrange so that 
u + (x) > u~(x), and v u {x) is pointing from the side of u~(x) to the one of u + (x). 

The measure Du can be canonically decomposed as the sum of three mutually 
orthogonal measures: the Lebesgue part Vu ■ C n , the jump part (u + — u~)v u ■ 
Ti n ~ 1 LSti, and a remainder, called Cantor part, which is singular but does not 
charge any W n_1 -finite set. 

The space SBV(£l) of special functions of bounded variation is given by all 
functions u 6 BV(£l) for which the Cantor part of the derivative vanishes, i.e., 



A subset E of Q has finite perimeter (in Q) if the distributional derivative DIe 
of its characteristic function is a bounded vector measure on f2. In this case, the 
measure theoretic boundary d*E is the singular set of 1e, while the inner normal 
v Q e is the associated normal vectorfield u 1e . Both the Lebesgue and the Cantor 
parts of DIe vanish, i.e., 



Du = Vu-C n + {u + - u~)v u ■ n n ~ l ^Su . 



(2.2) 



D1 E = v 



■H 



■n—l 



(2.3) 
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A general form of the divergence theorem 

When looking for calibrations for a given function, it is often convenient to 
consider also vectorfields which are not regular. In doing so, however, we face 
some technical difficulties. 

First of all, the first identity in (1.5) depends on the divergence theorem, and 
may not hold when <j) is divergence-free (in the sense of distribution) but not 
continuous, because the flux of such a vectorfield through a given surface is not 
well-defined. To solve this problem, we must assume a certain regularity in (f>. 

Definition 2.1. - We say that a vectorfield <j) on a subset E of W 1 is ap- 
proximately regular if it is bounded, and for every Lipschitz hypersurface M in 
R n there holds 

aplim [(f>(y) ■ v M (x)] = <f>(x) ■ v M {x) for U n ~ l -a.e. x € M n E , (2.4) 

where v M (x) denotes the (unit) normal to M at x. 

Remark 2.2. — If (j) is approximately regular, then (2.4) can be extended to 
every rectifiable set M, v M being now understood in the approximate sense. If <f> 
admits traces 4> + and 4>~ on the two sides of M (defined in the same way of the 
traces u + and u~), then (2.4) is equivalent to the compatibility condition 

</>• v M = </,+ ■ v M =</,-■ v M F-^ainMnfi, (2.5) 

which links the pointwise values of 4> on M with the values of the traces. 

In particular, if <j> is (approximately) continuous H n_1 -almost everywhere on 
E, then it is approximately regular. Similarly, if (j) is (approximately) continuous 
on the complement of a rectifiable set S, then <f> is approximately regular if and 
only if (2.4) holds for M := S. 

Remark 2.3. - If has the special form 4> := (0, ...,0, ip) where V = 
V>(xi, • • • ,x n ) is a bounded real function which is continuous in the variable x n , 
then (j) is approximately regular. Take indeed a Lipschitz surface M, and let Mq 
be the subset of all points x 6 M such that the n-th component of v M (x) vanishes. 
Then equality (2.4) obviously holds for all x 6 Mq. To prove that it also holds for 
7Y n_1 -a.e. point in M\Mq it suffices to notice that ip is approximately continuous 
in the complement of a set of type JVxR, where N is a £ n_ ^negligible subset of 
R n_1 . Thus 4> is approximately continuous at all points except those in N x R, 
which form an H n_1 -negligible subset of M \ M by the area formula (see [23], 
Theorem 3.2.22, or, for a slightly less general statement, [6], Theorem 2.91). 

We can now state a refined version of the classical divergence theorem (the 
proof is postponed to the appendix) . 
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LEMMA 2.4. - Let £1 be an open set in R™ with Lipschitz boundary, <J) an 
approximately regular vectorfield on Q, and u a function in BV(Q). Assume 
moreover that divcp € L°°(n) and ucp € L 1 (dQ,Ti n ~ 1 ). Then 



where u QVl is the inner unit normal to dVt, and in the last integral u stands for 
the trace of u on dQ. 

Notice that the condition ucft € L l (dVL,H n ~ l ) is always satisfied when dQ is 
bounded, because in this case the trace of u belongs to L l (dQ,H n ~ 1 ). 

Remark 2.5. - As shown in [7], for a Borel vectorfield (p with divergence in L 1 
it is possible to give a functional definition of the trace of the normal component 
of (f) on any Lipschitz surface M, and precisely as a continuous extension of the 
trace operator for regular vectorfields. When M is the boundary (of some set Q), 
this notion of trace automatically satisfies the divergence theorem (in the sense 
that formula (2.6) holds for every function u of class C*), but since modifying 
in a Lebesgue-negligible set affects neither this trace nor the distributional 
divergence, one can easily produce examples where the trace does not agree with 
the normal component of <fi on M. Thus Lemma 2.4 shows implicitly that this 
is never the case when <f> is approximately regular. One may wonder if every 
vectorfield with divergence in L°° agrees, up to a Lebesgue-negligible set, with 
an approximately regular one. Unfortunately the answer is negative, even for 
divergence-free vectorfields (cf. [7], example after Proposition 2.1). 

Going back to the first identity in (1.5), we remark that verifying that a 
vectorfield 4> is divergence-free is relatively easy when (ft is of class C 1 because 
the distributional divergence agrees with the classical one, which can be explicitly 
computed. If <f> is piecewise C 1 , the task is slightly more difficult, and can be 
carried out in many concrete cases (see Sections 4 and 5) with the help of the 
following lemma (the proof is postponed to the appendix). 

Lemma 2.6. - Let (f) be a bounded vectorfield on an open set 17 C R n , and 
assume that there exist a closed set S and a function f £ Lf oc (£l) such that 
divcfr = f in the sense of distributions onfl\S. Then the identity div0 = / holds 
also on Q if S can be written as S := So U Si, with So an H n ~ 1 -negligible closed 
set and Si a (possibly disconnected) Lipschitz hypersurface, and (f> satisfies (2.4) 
for M := S t and E := ft. 

Remark 2.7. — The point of this lemma is roughly the following: since the 
divergence is a first order differential operator, div^> cannot "charge" any set of 
codimension larger than 1, and therefore So can be safely removed. On the other 
hand, the part of divcp supported on the hypersurface M is given by the difference 
of the traces (whenever defined) of the normal components of cf> on the two sides 
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of M, which happens to vanish if (2.4) holds, and then we are allowed to neglect 
Si too. 

Remark 2.8. - Lemma 2.6 will be often applied in the following forms. 

(a) Suppose that is a bounded vectorfield on S7, continuous on Q \ (So U Si), 
and divergence-free on Q \ (So U Si), with So, and Si given as above. If (p satisfies 
(2.4) with M = Si, then cp is approximately regular on O and divergence-free on 
n (cf. Remark 2.2). 

(b) Let be given, for j = l,...,m, pairwise disjoint Lipschitz open sets Qj 
whose closures cover Q, and approximately regular, divergence-free vectorfields 
4>j on Qj. Let <p> be any vectorfield on Q which agrees at any point with one of 
the (pi (hence it is uniquely determined at least on the union of all f^). Then 
4> is approximately regular and divergence-free provided that the vectorfields <pi 
satisfy the compatibility conditions 

fii ' U an = $3 ' u dn H n ~ 1 -a.e. on dfli n dClj, 

which are equivalent to the compatibility condition (2.5) for (p. 

The complete graph of an SBV function 

We fix now some notation and state some results which are more specific to this 
paper. In the following O is a fixed bounded open subset of R n with Lipschitz 
boundary, and u m is its inner unit normal. The letter x usually denotes the 
variable in (or R n ), while t is the variable in R; U is an open subset of Q x R, 
with Lipschitz boundary whose closure can be written as 

U ■= {(x,t) G H x R : n(x) < t < r 2 (x)} , (2.7) 

where the functions ti,T2 : fi — > [— cx),+oo] satisfy n < ti- The letter <p> denotes 
a bounded vectorfield defined on (a subset of) R ra x R, with components cf> x G R™ 
and 4> l G R. Notice that div^ = d\v x 4) x + dt4> 1 , where div^ is the (distributional) 
divergence with respect to x and d t the (distributional) derivative with respect 
to t. 

Definition 2.9. — For every function u G BViVt), let l u be the characteristic 
function of the subgraph of u in Q x R, namely l u (x,t) := 1 for t < u(x) and 
l u (x,t) := fort > u(x). The complete graph ofu, denoted byTu, is the measure 
theoretic boundary of the subgraph ofu, i.e., the singular set of l u . 

Since the subgraph of u has finite perimeter in S7xR (see, e.g., [29], Proposition 
1.4), the definition of Tu is well-posed. Moreover Dl u = v Tu -H n LTw (cf. (2.3)), 
where v Vu is the inner unit normal of the subgraph of u. Therefore the flux 
through Fu of any vectorfield <f> on Q x R is given by the integration of (p with 
respect to the vector measure Dl u , that is 

/ <p-v ru dH n = [ 4>- Dl u . (2.8) 
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An alternative way to compute this flux is given by the following lemma (the 
proof is postponed to the appendix) . 



Lemma 2.10. - Let u be a function in SBV(Q) and let 
defined at least on Tu. Then 



be a vectorfield 



/> 

JTu 



u r dK 

1 u 



■ / y> x (x,u) ■Vu — (p t (x,u)]dx 
Jn 



+ 



Su 



(2.9) 



<j) x (x,t) dt 



■ u u dTC 1 - 1 



where u, u^, Vu, and v u are always computed at x. 



Remark 2.11. - Formula (2.9) corresponds to a decomposition of the deriva- 
tive of l u , or, better, to a decomposition of the complete graph Tu as union (up 
to 7Y ra -negligible sets) of a "regular" part - namely the set of all points (x,u(x)) 
such that u is approximately continuous at x, and has approximate gradient 
Vu(x) - and a "vertical" part - namely the set of all points (x, t) with x € Su 
and t 6 (u~ (x) , u + (x)) . Note that for a general BV function there would be an 
additional subset of Tu, corresponding to the Cantor part of Du. 

The following version of the divergence theorem (cf. Lemma 2.4) yields the 
first equality in (1.5) (the proof is postponed to the appendix). 

Lemma 2.12. - Let be given two functions u and v in BV{Q) whose com- 
plete graphs lie in U , and an approximately regular vectorfield <f> on U which is 
divergence- free in U. Then 



J ■ 



u r dn n 

1 u 



u r dU n 

1 v 



(2.10) 



: (x,t)dt 

(where, in the last integral, u and v stand for the traces on dCl). 



3. Calibrations for free discontinuity problems 

In this section we introduce a calibration principle for a wide class of free 
discontinuity problems, expanding the basic idea explained in the introduction. 
We begin with the case of the Mumford-Shah functional, with or without the 
lower order term, and then we consider more general functionals, possibly with 
discontinuous integrands, which include some interesting functionals considered 
in minimal partition problems. 

Throughout this section Q is a bounded open subset of R n with Lipschitz 
boundary, U is an open set with Lipschitz boundary contained in Q, x R which 
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satisfies (2.7); u always denotes a function in SBV(Cl). The functional F(u) and 
F (u) are given in (1.1) and (1.2), respectively, where Su and Vu are now defined 
as in Section 2, a > and j3 > are fixed constants, and g belongs to L°°(Cl). 
Note that the functional F is the particular case of F corresponding to a = 1 
and (3 = 0. In the following definition we fix some terminology about minimizers 
of F (which also applies to any other functional on SBV). 

Definition 3.1. - We say that a function u is an (absolute) minimizer of F if 
F{u) < F(v) for allv G SBV(Cl), while u is a Dirichlet minimizer if F(u) < F(v) 
for all v G SBV(Cl) with same trace on dCl as u (that is, with same boundary val- 
ues as u). We say that u is a U -minimizer if the complete graph of u is contained 
in U and F(u) < F(v) for all v G SBV (CI) with complete graph contained in U, 
while u is a U -Dirichlet minimizer if we add the requirement that the competing 
functions v have the same boundary values as u. 



Calibrations for the Mumford-Shah functional 



We generalize the idea described in the introduction and provide sufficient con- 
ditions for [/-minimality and {/-Dirichlet minimality with respect to the Mumford- 
Shah functional F (or Fo). We begin with the following key lemma. 

Lemma 3.2. - Let F be defined as in (1.1) for some a > 0, (3 > 0, and let (ft 
be a vectorfield on U which satisfies the following assumptions: 

(a) ^(x,t) > \\(j) x {x,t)\ 2 -[3(t-g) 2 for C n -a.e. x G Cl and every t G [ti,t 2 ], 

rto 



(b) 



(x,t)dt <afor7i n -a.e. x G Cl and every ti, t 2 G [ti, t 2 ], 



where the functions T\ and r 2 are defined by (2.7) and, like g, are computed at x. 
Then for every u such that Tu C U we have 



F(u)> [ 0-v ru dH r 



(3.1) 



Moreover, equality holds in (3.1) for a given u if and only if 

(a') (j) x (x,u) = 2Vu and ^(x,u) = \Vu\ 2 - (3(u - g) 2 for C n -a.e. x G Cl, 

(b') / (j) x (x,t)dt = au u for H 71 ' 1 -a.e. xG Su, 

J u~ 

where u, u* 1 , Vu, v u , and g are always computed at x. 

PROOF. Take u such that Tu C U. We recall that by Lemma 2.10 



/ 

JVu 



[4> x (x,u) ■ Vu - 0*(x,n)] dx 



+ 



Su 



(3.2) 



: (x,t) dt 



■ v^dH 11 - 1 . 
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It is an elementary fact that for every £, 77 6 R™ we have £ ■ 77 — ||£| 2 < |r/| 2 , and 
equality holds if and only if £ = 277. Hence, setting £ := (f) x (x,u) and 77 := V-u, 
and taking (a) into account, we obtain that, £ n -a.e. on Q, 

c/) x (x,u) ■ Vu-^(x,u) < (f> x (x,u) ■ Vu- ^\(p x (x,u)\ 2 + (3{u- g) 2 

< \Vu\ 2 + P(u - g) 2 

and consequently 

I [<j) x (x,u) ■ Vu- 0*(x,u)] dx < / [|Vu| 2 + /3(tj- 5 ) 2 ] dx . (3.3) 

Moreover, equality holds in (3.3) if and only if <j) x (x,u) = 2Vu and </>*(x,u) = 
||^> x (x,'u)| 2 - (3(u - g) 2 = \Vu\ 2 - P(u - g) 2 for a.e. x G ft, which is (a')- 

As for the second integral in the right-hand side of (3.2), condition (b) above 
implies 



"(x,t)dt 



■ Vu < 



"(x,t)dt 



< a W n-1 -a.e. on Su, 



and then 



(x,t) dt ■ v u dH 71 ' 1 < aH n - l {Su) . 



(3.4) 



Moreover it is clear that equality holds in (3.4) if and only if (b') is satisfied. 

Inequality (3.1) follows now from (3.2), (3.3), and (3.4), as well as the rest of 
the statement. □ 

Remark 3.3. - Let Q be the class of all vectorfields <fi on £1 x R, not necessarily 
bounded, which satisfy assumptions (a) and (b) of Lemma 3.2 with U := fl x R. 
It can be easily proved that for every u in SBV(£l) there exists <f> G Q which 
satisfies assumptions (a') and (b') for u, so that equality holds in (3.1). Starting 
from this, one can also find vectorfields <j) 6 Q of class C\ such that the value of 
the flux in the right-hand side of (3.1) is arbitrarily close to F(u). Taking (2.8) 
into account, we can thus prove that for all u G SBV{VL) there holds (cf. (1.6)) 



F(u) = sup 



sup 



Dl, 



(3.5) 



Moreover one can show that for any function u which is in BV(Q), but not in 
SBV(Cl), the last two terms in (3.5) are equal to +00. Since every integral of the 
form Jq xR ■ Dl u , with (j> of class C\ on VL x R, is continuous with respect to the 
weak* topology of BV(Q), formula (3.5) shows that the Mumford-Shah functional 
F, extended to +00 to the rest of BV(Q), is weak* lower semicontinuous. The 
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same argument can be easily applied to the general functionals considered in 
the next subsection, providing another proof of the well-known compactness and 
semicontinuity results in SBV due to L. Ambrosio (see [3], or [6], Sections 4.1 
and 5.4). 

Theorem 3.4. - Let u be a function with complete graph contained in U , 
and assume that there exists an approximately regular vectorfield <f> on U which is 
divergence-free on U and satisfies assumptions (a), (b), (a'), and (b') of Lemma 
3.2. Then u is a Dirichlet U-minimizer of F. If, in addition, the normal compo- 
nent of 4> at the boundary of U x R vanishes, namely 

^ . u dQ = H n -a.e. on (dfi x R) n dU , (3.6) 

then u is also an absolute U-minimizer of F. 

PROOF. Let v be a function in SBV (SI) such that v = u on dfl and Tv C U. 
Then 

F(v) > f 4>-v Vv dU n = [ <f>- u Tu dU n = F(u) . (3.7) 

JTv JFu 

Here, the first inequality and the last equality follow from Lemma 3.2, while the 
first equality follows from Lemma 2.12. We have thus proved that u is a Dirichlet 
[7-minimizer of F. Viceversa, assuming (3.6) we obtain that the first equality in 
(3.7) holds even if the traces of v and u on dSl differ, which proves that u is an 
absolute ^-minimizer of F. □ 

Definition 3.5. — We call the vectorfield 4> in the first part of Theorem 3.4 
a Dirichlet calibration for u on U (with respect to F). If <f> satisfies the additional 
assumption (3.6), then we call it an absolute calibration. 

When [/ := H x R we omit to write it. When it is clear from the context, we 
may also omit to specify the functional, the set U, and whether the calibration 
is Dirichlet or absolute, and simply say that is a calibration for u, or that <p 
calibrates u. 

Remark 3.6. - If <f> is an absolute calibration for u, then it is also an absolute 
calibration for every other minimizer. Indeed, if F(v) = F(u), the first inequality 
in (3.7) must be an equality, and by Lemma 3.2 this means that <p satisfies as- 
sumptions (a') and (b') for v too. Similarly, if is a Dirichlet calibration for u, 
then it is also a Dirichlet calibration for any other Dirichlet minimizer with the 
same boundary values as u. 

This fact can be sometimes used to prove that the minimizer is unique. 
For instance, if <fi calibrates a function u with a negligible singular set (i.e., 
TL n ~ 1 {Su) = 0), and the inequality in assumption (b) is always strict, then we 
deduce that assumption (b') can only be satisfied by functions with negligible 
singular sets, and therefore all minimizers should have this property. But on this 
class the functional F is strictly convex (for (3 > 0, and even for (3 = in case of 
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Dirichlet minimizers), and therefore the minimizer must be unique (see Remarks 
4.7, 4.11, 5.2, and Paragraphs 5.4 and 5.5). 

Remark 3.7. - The functional F in (1.2) is obtained by setting a := 1 and 
(3 := in the definition of F. In this specific case, assumptions (a), (b), (a'), and 
(b') in Lemma 3.2 become 

(a) 4> l {x,t) > j\(j) x (x,t)\ 2 for £ n -a.e. x G $7 and every t <E [ti,t 2 ], 

(b) / 4> x (x,t) dt < 1 for H n ~ 1 -a.e. x G and every ti,t 2 G [ti,t 2 ], 
(a') 0*(a;,u) = 2Vu and ^{x,u) = |Vu| 2 for £ n -a.e. xe!l, 

(b') / x (x, t) dt = v u for H n_1 -a.e. x G Su. 

Remark 3.8. — It must be noticed that, given a boundary value w in 
L 1 (d£l,H n ~ 1 ), the Dirichlet problem 

min {f (u) : u G SW(fi), u = w r^-a.e. on 0ft, Tu C I/j (3.8) 

may not have a solution, even for a very regular to, due to a lack of continuity 
of the trace operator on SBV(Q). Therefore problem (3.8) is usually replaced by 
the relaxed problem 

in |f (u) + H n ~ x ({x € dQ : u{x) ^ w(x)}) : u G SBV(Q), Tu C I/} (3.9) 



mm 



(where in the second term u denotes, as usual, the trace of u on <9f2). A variant 
of the standard lower semicontinuity and compactness theorems in SBV(Q) (see 
[14] or [8]) shows that problem (3.9) has always a solution. 

The calibration method applies to problem (3.9), too. In this case calibrations 
are approximately regular vectorfields ^ on J] x R which are divergence-free, 
and satisfy conditions (a), (b), (a'), (b') of Remark 3.7 and, in addition, the 
following two conditions (which may be viewed as extensions of (b) and (b') to 
the boundary): 

(c) / 4> x (x,t) dt < 1 for H"- _1 -a.e. x G d£l and every s G [ti,t 2 ], 

J W 

(c') / (f) x (x, t) dt = u 9Q for H n_1 -a.e. x 6 dQ with u(x) + w(x), 

J w 

where u and w are computed at x, as well as t\ and r 2 . Indeed, arguing as in 
Lemma 3.2 we obtain 

Fo^+H 71 ' 1 ({x G dn : v(x) + w(x)}) > 

>[ 4>-v Tv dH n + [ \f (f> x (x,t)dt\-u dn dn n - 1 ( "' l0) 

Jrv Jon l Jw J 

for every v G SBV(Q), with equality for v = u, moreover the right hand side of 
(3.10) does not depend on v (apply Lemma 2.12 taking into account that <p IS 
divergence free). 
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Calibrations for general functionals 

The method of calibrations can be easily adapted to a larger class of functionals 
on SBV(Q). Take indeed 

f f(x,u,Vu)dx + [ i\){x,u- ,v u )dH n - x , (3.11) 
Jn J Su 

where / : fixRx R n — > [0, +oo] andi|>:fixRxRx S™" 1 -» [0, +oo]. We refer 
to [3] for general conditions on / and ip which imply the lower semicontinuity of 
the functional (3.11) and guarantee the existence of minimizers. However, lower 
semicontinuity is irrelevant for the calibration method. 

Let /* and d$f denote the convex conjugate and the subdifferential of / with 
respect to the last variable. We recall that the subdifferential of a function g : 
R n — > [0, +oo] at the point £ G R n is defined as the set of vectors r\ G R n such 
that g(£) + 1] ■ {( — < g{() for every ( G R n . It is well known that for every 
£, T) G R n we have the inequality £ • 77 — <?*(f?) < <?(£), and that equality holds if 
and only if 77 G (?£<?(£) . Using these properties we obtain the following variant of 
Lemma 3.2, whose proof is omitted. 

Lemma 3.9. - Let <j) be a vectorfield on U which satisfies the following as- 
sumptions: 

(a) t) > f*(x, t, 4> x (x, t)) for C n -a.e. and every t G [n, r 2 ] ; 

(b) / x (x, t) dt -v < ip(x, ti,t 2 , v) for Ti 11-1 -a.e. every v G S n_1 ; 

and every t\ < t 2 in [ti,T2]. 
T/ien /or every u wi/i complete graph contained in U we have 

> / ^ • ( 3 - 12 ) 

Moreover, equality holds in (3.12) /or a given u if and only if 

(a') (f) x (x,u) G d^f(x,u,Vu) and 0*(a;,it) = f*(x,u,(p x (x,u)) for C n -a.e. 

(b') / <f> x (x,t)dt ■i/ u = ip(x,u~,u + ,v u )forTC n ~ 1 -a.e.x£Su, 

J u~ 

where u, , Vu, and v u are always computed at x. 

Proceeding as in the previous subsection, one can prove the following analogue 
of Theorem 3.4. 

Theorem 3.10. - Let u be a function with graph contained in U. Assume 
that there exists an approximately regular vectorfield (j> on U which is divergence- 
free and satisfies assumptions (a), (b), (a'), and (b') of Lemma 3.9. Then u 
is a Dirichlet U-minimizer of ^> . If, in addition, the normal component of 4> 
on the boundary of Q x R vanishes, i.e., (3.6) holds, then u is also an absolute 
U -minimizer of ^> . 
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Calibrations for minimal partitions 

Besides F and Fq, interesting examples of functionals of the form (3.11) arise 
in different variants of the minimal partition problem. Let us consider the case 
where the number m of the elements of the partition is prescribed. To formulate 
the problem, we fix real numbers a\ < ■ ■ ■ < a m in an arbitrary way, and consider 
only functions u in the class BV(Cl, {ai}) of all u 6 BV(Cl) which take only 
these prescribed values. The corresponding level sets Ai := {u = o^}, sometimes 
called phases, form a partition Q; for % < j, the set SV,- of all x 6 Su such that 
u~{x) = cii and u + {x) = aj is called the interface between the phases Ai and A,-, 
and is oriented by the normal i/jj pointing from Ai to Aj (hence = v u ). 

We consider functionals of the form 



where ipij : ft x S ra_1 — > [0, +oo]. Notice that the weights tpij may depend on the 
point at which two phases meet and on the normal to the interface at that point. 
A partition (A\, . . . , A m ) is said to be a Dirichlet minimizer of T if it minimizes 
T among all partitions (Bi,...,B m ) such that, for every i, the characteristic 
functions of Ai and Bi have the same trace on dVi. Notice that all these notions 
do not depend on the particular choice of the numbers a\, . . . , a m . 
We now define a functional ^ of type (3.11) by setting 



Since every function u in BV(Q,{ai}) belongs to SBV(Q), and V-u = £ ra -a.e. 
in 17, one easily checks that ^(u) is finite only if u belongs to BV(Q, {ctj}), and in 
this case ^(u) = T(Ai, . . . , A m ), where (A\, . . . , A m ) is the partition associated 
with u. Hence a partition (A\, . . . , A m ) is a Dirichlet minimizer for T if and only 
if the corresponding function u is a Dirichlet minimizer for fy, which according to 
Theorem 3.10 is implied by the existence of a calibration. 

For this particular choice of a calibration for u 6 BV(Q,{ai}) is an ap- 
proximately regular vectorfield 4> on VL x R which is divergence-free and satisfies 
the following properties (cf. Lemma 3.9): 

(a) 0*(x, aj) > for £ n -a.e. x 6 Q and every i, 



J r (A 1 ,...,A rn ) = ^2 '.';/!•''. I', j ) (III 



(3.13) 



i<3 ij 



f(x,t,0 : = 
ip(x,t 1 ,t 2 ,u) : = 



if t € {ai} for some i and £ = 0, 
-f-oo otherwise, 

ipij(x, v) if t\ = ai and t 2 = aj for some i < j, 
+oo otherwise. 




-a.e. x £ Q and every v G S 



n— 1 



and i < j, 

(a') 0*(x, aj) = for £ n -a.e. x G A, and every i, 
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(b') / cf) x (x,t) dt ■ Uij = i/)ij(x, Uij) H n_1 -a.e. on Sij for every i < j. 

Ja,i 

These calibrations can be re-written in a different and more interesting way: 
let us set, for i = 1, . . . , m, x € fi, 

4>i{x) := \ (j) x (x,s) ds . 



ai 



One can check that the vectorfields are approximately regular on VL and have 
divergence in L°°(Q) (more precisely, div0j(x) = (f) x (x,ai) — (f) x (x,ai) - apply, 
e.g., formula (2.6) with f2 replaced by ft x (ai,aj), and u any smooth function on 
x (01,0^) which depends only on x) and satisfy the following properties: 

(c) div^i > div<^ for £ n -a.e. x £ Ai for every j ^ i, 

(d) (<f)j — 4>i) ■ v < ipij(x, v) for H n ~ 1 -a.e. x 6 Q and every v 6 S n_1 , i < j, 
(d') (4>j — 4>i) ■ uij = i/)ij(x, Uij) for H n_1 -a.e. x € S'jj and every i < j. 

Conversely, given approximately regular vectorfields 4>i ° n ^ with divergence in 
L°°(Q) which satisfy (c), (d), (d'), we define a vectorfield cf> on Q x R as follows: 
we take smooth non-negative functions cr, with support included in (cij, ctj+i) and 
integral equal to 1, and set 

4> x {x,t) := (jj(t) (0 i+ i(x) - 4>i{x)) iov x e CI, ^ < t < a i+1 , 

then we take 0* so that 

(/>*(#, := for x € Aj, 

9 t 0'(x,t) := <n(t) (div^(x) - div0 l+1 (x)) (3-14) 

for x G 0, < i < aj+i; 

the definition is completed by setting <p(x,t) := (f){x,a\) for i < ai and c/)(x,t) := 
(j)(x,a m ) for t > a m . 

Then is divergence free in Q x R by construction, and one can easily check 
that it satisfies assumptions (a), (b), (a'), and (b') above. Moreover, since each fa 
is approximately regular on O, one can verify that (cf> x , 0) is approximately regular 
on Q x R, and the same holds for (0, (j) 1 ) by Remark 2.3. Hence <j> is approximately 
regular too. 

We have thus proved the following result. 

Theorem 3.11. - Let (Ai, . . . ,A m ) be a partition of Q, and assume that 
there exist approximately regular vectorfields fa , . . . , (p m on f2 with divergences in 
L°°(fi) which satisfy assumptions (c), (d), and (d') above. Then (Ai, . . . ,A m ) is 
a Dirichlet minimizer of the functional T in (3.13). 

Remark 3.12. — A particularly relevant example of functional of type (3.13) 
is the "interface size", which is obtained by taking ipij = 1 for all % < j. In this 
case assumptions (d) and (d') above reduce to 
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(d) \4>j{x) — <f>i(x)\ < 1 for every x £ and every i < j, 

(d') (f)j(x) — <f>i(x) = Vij{x) for H n_1 -a.e. x € and every i < j. 

Calibrations of this type have already been introduced in [34] and [11] 
as "paired calibrations". More precisely, a paired calibration for a parti- 
tion (Ai, . . . , A m ) is an ordered m-uple of approximately regular vectorfields 
(j>i, ... , 4> m on f2 which are divergence-free and satisfy assumptions (d) and (d') 
above. We notice that the assumption that the vectorfields 4>i are divergence-free 
is stronger than (c), and indeed our definition allows in principle for a larger class 
of calibrations. 

Among other applications, in [34] it is shown that in any dimension n the 
partition of a regular simplex in R n given by the n + 1 simplices spanned by one 
face and the centre is a Dirichlet minimizer of the interface size (and the paired 
calibration consists simply of n + 1 constant vectorfields which are orthogonal 
to the corresponding faces). For n = 3 this statement was first shown in [38] 
with a (relatively) long proof. In [11] it is shown that, unlike what happens in 
dimension 3, the partition of a hypercube in R n , n > 4, given by the 2n simplices 
spanned by one face and the centre is a Dirichlet minimizer of the interface size. 
In both papers the theory is also extended to cover more general functionals. For 
further examples and results, see also [13], [12]; other references are included in 
the survey [32]. 

4. Applications to the homogeneous 
Mumford-Shah functional 

In this section we give some examples of Dirichlet minimizers of the homoge- 
neous Mumford-Shah functional F . We begin with a few remarks which may be 
useful when constructing calibrations. 

Remark 4.1. - By a simple truncation argument, to prove that a function 
u : f2 — > [rn, M] is a (Dirichlet) minimizer for Fq it suffices to show that Fq(u) < 
Fq(v) for all competitors v such that m < v < M. Thus it is enough to show that 
u is a Dirichlet {/-minimizer, with U := ft x (m, M). In the following we often 
tacitly assume this principle, and construct calibrations infix [m, M] instead of 
x R. Notice, however, that a calibration <j> on Cl x [rn, M] can be extended to 
x R in a rather simple way: it suffices to set (fi(x,t) := (0, 0* (x, mj) for t < m 
and <l>{x,t) := {0,^{x,M)) for t > M (cf. Remarks 2.3 and 2.8(b)). 

The same conclusion holds for F if g satisfies m < g < M too (but may fail 
for a functional of the general form (3.11), due to lack of suitable truncations). 

Remark 4.2. — We can construct divergence-free vectorfields on an open set 
A C Q x R using fibrations of A by graphs of harmonic functions. More precisely, 
given harmonic functions u\ whose graphs are pairwise disjoint and cover A, for 
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all (x, t) 6 A we set 

4>(x,t) := (2Vu A (x),|Vu A (x)| 2 ) , (4.1) 

where A = A(x,t) is taken so that t = u A (x). Thus cj> is a vector-field on ^4 which, 
by construction, satisfies assumption (a) of Remark 3.7, and assumption (a') for 
every u\. 

We prove that 4> is divergence-free under the additional assumption that the 
function u(x, A) := u\(x) is of class C 1 and d\u(x, A) ^ for every (x, A), which 
implies that the parameter A can be (locally) chosen so that it depends on x and 
t in a C 1 fashion. Then we get 

div</> = 2A x u + 2d\V x u ■ V X X + 2V x u ■ d\V x ud t X 

= 2d x V x u ■ (V X A + V x u d t X) . (4.2) 

On the other hand, deriving the identity u(x,A(x,t)) = t with respect to x and 
t we get V x u + d\uV x X = and d\ud t \ = 1, respectively. This implies that the 
last factor in (4.2) vanishes, and thus 4> is divergence-free (to make this argument 
work, we need that V x u is of class C 1 in (x, A), which can be derived by the fact 
that each function u\ is harmonic). 

In Paragraphs 4.6 and 4.10 below we apply this idea by embedding a harmonic 
function that we intend to calibrate into a family of harmonic functions whose 
graphs fibrate A := Q x (m,M), and taking cp as in (4.2). To show that is a 
calibration we will have only to verify assumption (b) of Remark 3.7. 

Remark 4.3. — The construction described in the previous remark is a par- 
ticular case of a classical result about extremal fields of scalar functionals (see, 
e.g., [1], Section 4). Let be given an open subset A of £1 x R which is covered by a 
family of pairwise disjoint graphs of solutions u\ of the Euler-Lagrange equation 
associated with the functional f n f(x, u, Vu) dx, that is 

div(<9 e /(x,n, Vu)) = d t f(x,u,Vu) . 

Here t and £ denote the second and third variable in the argument of f(x,t,£), 
corresponding to u and Vu, and we assume that / is of class C 2 in (i,£) and 
convex in £. For all (x, t) € A let 

(f> x (x,t) := d i f(x,u x (x),Vu x (x)) , 

(f>*(x,t) := dzf(x,u\(x),Vu\(x)) ■ Vn A (x) - /(x, u x (x), Vn A (x)) 

where A = X(x,t) is taken so that t = u\(x). Then </> satisfies by construction as- 
sumption (a) of Lemma 3.9, and assumption (a') for every function u\. Moreover 
one can prove that <p is divergence-free under the additional assumption that the 
function u(x, A) := u\(x) is of class C 2 in (x, A) and d\u ^ (cf. [1], Theorem 
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4.6). If, in addition, <j) satisfies assumption (b) of Lemma 3.9, then it is a cali- 
bration for each u\ relative to the functional ^ in (3.11), and hence each u\ is a 
[/-Dirichlet minimizer of Vl/ provided that U is contained in A. 

Remark 4.4. - The construction described in the previous remarks cannot 
be really used for absolute minimizers, i.e., when no Dirichlet boundary condition 
is imposed. Indeed, calibrations for absolute minimizers should have vanishing 
normal component at the boundary of x R (see Theorem 3.10), and this holds 
true for the vectorfield (f> in (4.3) if and only if the functions u\ satisfy the natural 
boundary condition d^f(x,u\(x),Vu\(x)) ■ v dn (x) = for x £ <9H But this 
means that all u\ solve the Neumann problem associated with the functional 
J Q f(x, u, Vu) dx, while in general we cannot expect a one-parameter family of 
solutions for such a problem. For instance, if this functional is strictly convex 
(and this is indeed the case for the regular part of the Mumford-Shah functional F 
when (3 > 0), then the associated Neumann problem admits at most one solution. 

Remark 4.5. - For n = 1 and Q = (a, b), the equation div^ = on x R, 
coupled with the identity 0* = j((f) x ) 2 , reduces to the first order equation 

d x (p x + \<p x d t <t) x = . (4.4) 

It easily follows from the method of characteristics that 4> is a C 1 solution of 
(4.4) in Q x R if and only if every level set {<jf = s} is composed of straight 
lines with slope s/2 (intersected with $7 x R). In other words, for n = 1 all C 1 
divergence-free vectorfields 4> on x R which satisfy = \ {(f> x ) 2 (cf. conditions 
(a) and (a') in Remark 3.7) are associated with a fibration of O x R with graphs 
of affine -i.e., harmonic- functions as in Remark 4.2. 

For the rest of this section, calibrations are always intended as Dirichlet cal- 
ibrations for F , in the sense of Remark 3.7. We begin with a discussion of 
some one-dimensional examples. Of course, in these examples minimality can be 
easily checked by direct computations, and there would be no need for calibra- 
tions. Nevertheless, it is instructive to see what happens, and moreover some 
one-dimensional constructions are carried over to higher dimensions (cf. Para- 
graphs 4.10 and 4.12). 

4.6. Affine function in one dimension. - Let Q be the open interval (0, a) 
and let u be the linear function u(x) := Ax, with A > 0. It is easy to see that u 
is a Dirichlet minimizer of Fq if and only if 

a\ 2 < 1 . (4.5) 

In this case a calibration is given by the piecewise constant vectorfield: 

f(2A,A 2 ) if \x < t < kx + a), 
<Kx,t):=\) \ ) 2 " " 21 (4.6) 
^ (0, 0) otherwise. 
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Thus 4> satisfies assumptions (a) and (a') of Remark 3.7, and vanishes outside a 
stripe of constant height (in grey in Figure 2 below, on the left) which is arranged 
so that (b) holds and div</> vanishes (cf. Remark 2.8(b)), while (b') is trivially 
satisfied. 




-+-x 



Figure 2 



Another calibration is obtained by fibrating the rectangle U = (0, a) x (0, Aa) with 
affine functions as shown in Figure 2, on the right, and applying the construction 
of Remark 4.2: 



<f>(x,t) :-- 




if < t < Ax, 
if Ax < t < Aa. 



(4.7) 



It remains to check that assumption (b) is satisfied, which happens if and only if 
(4.5) holds. 

Remark 4.7. - If a\ 2 < 1, then both calibrations described in the previous 
paragraph satisfy the strict inequality in assumption (b) of Remark 3.7, i.e., 



(x, t) dt < 1 for every x G [0, a] and every t\, £2 £ R- 



By Remark 3.6, this shows that the function u(x) := Ax is the unique Dirichlet 
minimizer of F with u(0) = and u(a) = Aa. 

4.8. Step function in one DIMENSION. - In Paragraph 4.6, in the limit case 
aA 2 = 1 the linear function u(x) = Ax and any step function of the form u(x) := 
for < x < c and u(x) := Aa = ^/a for c < x < a (with < c < a) are both 
Dirichlet minimizers with the same boundary values. Hence both vectorfields 
(4.6) and (4.7) calibrate these step functions when A := l/yfa (cf. Remark 3.6). 
Furthermore, it is easy to check that they also calibrate any step function u given 
by u(x) := for < x < c, and u(x) := h for c < x < a with h > y/a. 

Remark 4.9. - When aA 2 > 1 the linear function u(x) := Ax is not a Dirich- 
let minimizer of Fq (a step function is preferable), but it is still a Dirichlet U- 
minimizer, when U is the stripe of all points (x, t) between the graph of Ax — 
and Ax + ^ . A calibration is given by <p(x, t) := (2A, A 2 ). 
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^-neighborhood of Zh 




-heigth=Ve ! u 
slope=(4Ve) _1 



Figure 3 

Conversely, when h < y/a, the step function u in Paragraph 4.8 is no longer a 
Dirichlet minimizer, but it is Dirichlet [/-minimizer when U is an e-neighbourhood 
of the complete graph of u (in grey in Figure 3) and e satisfies \\fs + 2e < h. A 
calibration is given by the piecewise constant vectorfield which vanishes outside 
the white parallelogram in Figure 3, and is equal to ^) inside. 

4.10. Harmonic functions in dimension n. - Let u be a harmonic function 
on f2. Since u is a Dirichlet minimizer of |Vu| 2 , it is natural to ask when it is 
also a Dirichlet minimizer of F . As pointed out by A. Chambolle, this happens 
when 



osc u • sup |Vu| < 1 , 



(4.8) 



where osc u is the oscillation of u, namely the difference between the supremum 
M and infimum m of u (over In the one-dimensional case n = 1 this condition 
reduces to (4.5). 

A calibration can be constructed by analogy with (4.6); see Figure 4, on the 
left: 



f (2Vu(x), \Vu{x)\ 2 ) if \{u{x) +m)<t< \{u{x) + M), 
(f)(x,t) := < 

[ (0,0) otherwise. 



(4.9) 



Another calibration can be obtained, as the one in (4.7), by embedding u in a 
family of harmonic functions whose graphs fibrate the cylinder !lx[m, M]. More 
precisely we take the functions m + \(u — m) and M + \{u — M) with A ranging in 
[0, 1] (see Figure 4, on the right), and then the construction of Remark 4.2 gives 



(x,t) :-- 



! ^)^V U (x), {^_f\Vu(x)\ 2 ) if m < * < u(x), 



(2 



M-t 
M-u(x) 



V«(s),( 



M-t • 
M-u(x) • 



l \Vu(x)\ 2 ) \iu{x)<t<M. 



(4.10) 



One easily checks that both vectorfields are divergence-free (see Remarks 2.8(b) 
and 4.2), and that assumptions (a) and (a') of Remark 3.7 are satisfied, assump- 
tion (b') is always trivially satisfied, while assumption (b) holds if an only if (4.8) 
holds. 
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When (4.8) is not satisfied, u is still a Dirichlet {7-minimizer of Fq, where 
U is the set of all points (x,t) € fl x R which lie between the graph of 
u{x) — (4|Vn(x)|) _1 and u(x) + (4| Vu(x)|) _1 ; a calibration is given by <j)(x,t) := 
(2Vu(x),\Vu(x)\ 2 ). 

Remark 4.11. - If inequality (4.8) holds and u is not affine, then the maxi- 
mum principle implies that osc it - |Vit(x)| < 1 for every iGll, and therefore both 
calibrations constructed in the previous paragraph satisfy the strict inequality in 
assumption (b) of Remark 3.7. By Remark 3.6, this proves that the harmonic 
function u is the only Dirichlet minimizer of F with the same boundary values 
as u (cf. Remark 4.7). 

4.12. Step function in dimension n. - Let SI be a product of the form 
(0,a) x V, where V is a regular domain in R n_1 , n > 2, and let u be the step 
function given by u(x) := for < X\ < c, and u(x) := h for c < X\ < a, where 
X\ denotes the first coordinate of x, c G (0, a), and h > y/a. 

Using the results in Paragraph 4.8 and a symmetrization argument, it is easy to 
see that u is a Dirichlet minimizer of Fq. Calibrations can be constructed starting 
from the one-dimensional ones described in Paragraph 4.8. More precisely, we 
take the vectorfield on Q x R which is parallel to the {x\ , t)-plane and is given by 
formula (4.6) (or even (4.7)) with x replaced by X\ and A := \ j\fa. 

Remark 4.13. - The previous result can be restated by saying that a step 
function u with jump of height h along a hyperplane orthogonal to the direction 
e is a Dirichlet minimizer whenever h > \/a, where a is the diameter of the 
projection of £1 along the e axis. As in dimension one, any step function is a 
Dirichlet [/-minimizer, when U is an e-neighbourhood of the complete graph of u 
and e satisfies \\J~s + 2e < h (cf. Remark 4.9). However, unlike what happens in 
dimension one, u may be a Dirichlet minimizer of Fq even when the assumption 
h > \fa is not satisfied (cf. Paragraph 4.19 below). 

4.14. Triple junction in the plane. - Let Q := B(0,r) be the open 
disk in the plane with radius r and centred at the origin, and let (A 1: A 2 , A 3 ) 
be the partition of £1 defined as follows: A4 is the set of all x 6 Q of the form 
x = (pcos9, psinO), with |7r(i — 1) < 6 < Finally define u := a, in each Ai, 
where a\ < a 2 < as are distinct constants. 

Thus the singular set of u is given by three line segments Si^, 52,3, an d 
5*3,1 meeting at the origin with equal angles (see Figure 5, on the left), and it is 



20 



G. Alberti, G. Bouchitte, G. Dal Maso 



well-known that this is a minimal network, in the sense that the corresponding 
partition (A\, A 2 , A3) is a Dirichlet minimizer of the "interface size" functional 
(see Remark 3.12). Therefore it is natural to conjecture that, when the values of 
the constant ctj are sufficiently far apart, u is a Dirichlet minimizer of F too, that 
is, there is no convenience in removing part of the jump and taking a function 
with non- vanishing gradient. 

We prove this conjecture by calibration. Inspired by the constructions de- 
scribed in Paragraphs 4.8 and 4.12, we take e± := (±-\/3/2, — 1/2), A > 0, and 
set 

if \t - |(ai + a 2 ) - I x ■ e_| < ^, 
if \t-\(a 2 +a z )- fx-e+| < ^, (4.11) 
otherwise. 

Thus is piecewise constant, satisfies assumption (a) of Remark 3.7 by construc- 
tion, and vanishes out of two slabs of constant height ^ (see Figure 5, on the 
right). These slabs have been arranged in order to fulfill the following require- 
ments: 

(i) one slab is contained in Q x [en, a 2 ] and the other one in Q x [02, 0,3], so 
that assumption (a') of Remark 3.7 is satisfied; it is possible to construct 
such slabs if we can choose A so that dj + i 



(x,t) :-- 



(2Ae_,A 2 ) 
(2Ae+,A 2 ) 
(0,0) 



> Ar + t^, that is, if 



(4.12) 



(ii) the compatibility condition (2.5) is satisfied on the boundary of the 
slabs, so that 4> is approximately regular and divergence-free (cf. Remark 
2.8(b)); 

(iii) assumption (b') is satisfied for all points x in and 6*2,3, where v u 
coincides with e_ and e + respectively. 




first slab 
$=(2Ae + ,A 2 ) 



Q=B(0,r) 



second slab 
<H2Ae_,A 2 ) 



Figure 5 




t-a 



t-a. 



Moreover (b') holds also for x in S^i, because e_ + e + = v u . One also checks 
that the integral 4> x (x,t) dt can be always written as a linear combination 
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/i_e_ + /U_|_e_|_ with p± in [0, 1] (depending on x,ti,t 2 ), and since e + and e_ span 
an angle equal to 27r/3, this implies that the integral has modulus not larger 
than 1. Thus (b) holds, too. 

Remark 4.15. - When a 2 — a\ (or 03 — 02) is sufficiently small, u is not a 
Dirichlet minimizer. More precisely, if 02 - ox < (cf. (4.12)), a comparison 
function v with the same boundary values as u and such that F(v) < F(u) is 
given, in polar coordinates, by 

'5(01+02) if < 6 < |vr and p < r - d, 

2(01 + 02) + 5(01 — 02) \{P — r + d) if < 9 < g7r and p > r — d, 

v := < 

2(01 + a 2 ) + \{a 2 — ai) i(p — r + d) if |7r < 6* < g7r and p > r — d, 

k a 3 if |?r < 6> < 2vr, 

where d := (02 — a\)y/r (we leave the computations to the reader). 

4.16. Minimal partitions in dimension n. - One can generalize the exam- 
ple of the triple junction, and conjecture the following: if a partition (Ai, . . . , A m ) 
of is a Dirichlet minimizer of the "interface size" (see Remark 3.12) and u is a 
function which takes a constant value Oj on each Ai (with a\ < a 2 < ■ ■ ■ < a m ), 
then u is a Dirichlet minimizer of Fq when the values Oj are sufficiently far apart 
from each other. Unfortunately we can only prove this statement under two ad- 
ditional assumptions: 

(i) the partition (Ai, . . . , A m ) is not only minimal, but admits a paired cal- 
ibration in the sense of [34] and [11], namely there exist approximately 
regular, divergence-free vectorfields fa , . . . , <\> m on Q which satisfy as- 
sumptions (d) and (d') in Remark 3.12; 

(ii) for i = 1, ... , m — 1 there exist Lipschitz functions tpi : Q — > R which 
satisfy almost everywhere the first order equation 

• (&+i - <tn) = \\<t> i+ i - H 2 ■ (4-13) 

Adding, if needed, a constant to ipi, we may also assume that 

osc^i = 2||Vi||in/iy . (4.14) 

For i = 1, . . . , m — 1 we take slabs E/j, included in Q x (aj, Oj+i), of the form 

Ui = |(ar,t) : |t - |(oi + a m ) - A^^)! < 4^7} , (4.15) 
where the constants Ai will be chosen below. Then we set (cf. (4.11)) 

r {2\ i (4> i+1 (x) - Ux)), a? \ct> i+1 (x) - M*)\ 2 ) 

4>{x, t) := < if (x, t) G Ui for some i, 

I (0, 0) otherwise. 
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Taking into account assumption (d') in Remark 3.12 and the definition of the 
slabs Ui, one can easily check that assumptions (a), (a'), and (b') of Remark 3.7 
are satisfied. 

Let us check assumption (b). Taken t\ G [aj,ctj_|_i] and t<i G [ay, aj+i] for some 
the integral J t * 2 (p x (x,t)dt can be decomposed as the sum of the integrals 
on the (oriented) intervals [ii,ai+i], [a i+1 ,aj], and [aj,t 2 ], and hence it can be 
written as 

m((/> i+1 (x) - 4>i(x)) + (4>j(x) - 4> i+ i(x)) + fi 2 (4>j+i(x) - 4>j(x)) 

for suitable Hi,fJ,2 G [0,1]- But this sum can be reorganized as the difference 
between fj, 2 (pj+i(x) + (1 — (j, 2 )(pj(x) an d (J-i<fri(x) + (1 — fii)(pi + i(x). Therefore 
its modulus is the distance between two points in the convex hull of the vectors 
01 (x) , . . . , 4> m (x) , which has diameter 1 because of assumption (d) in Remark 
3.12, and (b) is proved. 

Since the vectorfields (pi are divergence-free and approximately regular by 
assumption, <fi is divergence-free and approximately regular within each slab (the 
approximate regularity of (cp x , 0) is immediate, that of (0, </>*) follows from Remark 
2.3), as well as in the interior of the complement of the union of all slabs. Thus 
<j) is divergence-free and approximately regular in Q x R if (and only if) the 
compatibility condition (2.5) is satisfied on the boundary of each slab (cf. Remark 
2.8(b)), which reduces to equation (4.13). 

Therefore we have constructed a calibration for u, provided that we can choose 
Aj so that the slabs U are contained in Q x (aj,aj + i), that is, 

— CLi .. .. 1 Xi 1 

and this can be done as long as 

a i+ i — a,i > \J 2 osc tpi for i = 1, . . . , m — 1. (4-16) 

Remark 4.17. - A paired calibration for the partition (Ai, A 2 , A3) de- 
scribed in Paragraph 4.14 is given by the constant vectorfields (f>\ := ^(v / 3, 3), 
(f>2 '■= |(— 2\/3, 0), 03 := ^(v / 3, — 3), and the linear functions tpi and ^2 with 
gradients \(— V3, — 1) and \(V3, — 1) satisfy equation (4.13); if we thus apply 
the construction of the previous paragraph, we obtain exactly the calibration 
described in Paragraph 4.14. 

Remark 4.18. - The first order equation (4.13) does not always admit so- 
lutions. For instance, since the derivative of tpi along the integral curves of 
the vectorfield </>j+i — (pi (i.e., the maximal solutions of the differential equation 
7 = <pi + i(j) — (pi('*/)) is always positive, when there exists a nontrivial closed inte- 
gral curve within Q, (4.13) admits no solution. On the other hand, if (pi + \ — (pi is 
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C 1 and nowhere vanishing, and all integral curves start and end at the boundary 
of and intersect a fixed (n — 1) -dimensional closed manifold M in ft which is 
transversal to the vectorfield fa + i — fa, i.e., a cross-section of the associated flow, 
then the method of characteristics provides a solution tpi to (4.13) of class C 1 . 

However, such a strong requirement on fa+i — (pi is far from being necessary. 
Not only there may exist Lipschitz functions fa which satisfy (4.13) almost every- 
where even if fa+\ — fa vanishes somewhere, but for the purposes of the previous 
construction we can even allow fa to be discontinuous along some integral curve 
7: in this case the boundary of the slab Ui in (4.15) is not just the union of the 
graphs of Xiipi + +02) + 4T7 and + + 0-2) — jtt, but there is an 
additional vertical piece contained in 7 x R. Yet the compatibility condition (2.5) 
is satisfied there, and then cj> is still divergence-free and approximately regular (cf. 
Remark 2.8(b)). 

4.19. More on the step function in the plane. - Let Q be the rectangle 
(—a, a) x (—b, b) in the plane, and let (Ai,A 2 ) be the partition of Cl given by the 
sets of all points x = (#1,2:2) such that x\ < and x\ > respectively. This 
partition is obviously minimal, and a paired calibration is fa := (0, 0), fa '■= (1, 0). 
Setting p + := (0,6) and p- := (0,-6), another paired calibration is given by 
fa := (0,0) and 



where p±,6± are the polar coordinates around the points p±; see Figure 6. 



In this function ipi which satisfies (4.13) almost everywhere is given by 




(— sin 6 + , cos 0+) for x € B(p+, b), 
(sin#_, — cos6L) for x € B(p_, b), 
(0, 0) otherwise, 



(4.17) 




Figure 6 




i(# + +7r/2)p + for x e B(p+,b), 
±(0_-7r/2)p_ for x G B{p_,b), 



otherwise, 
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and the construction in Paragraph 4.16, performed with some care because of 
the discontinuity of ip\ along the circles dB(p±,b), yields a calibration for the 
step function u which takes the value a\ on A\ and ai on A^-, provided that 
0-2 — o-i > VttB- Note that the calibration obtained in this way is defined on the 
whole stripe R x (—b, b) and does not depend on a. This extends the minimality 
result proved in Paragraph 4.12. 

4.20. More on the triple junction. - Let us apply again the con- 
struction described in Paragraph 4.16 to the situation described in Paragraph 
4.14, with Vt replaced by an e-neighbourhood of Su within the ball B(0,r) (in 
grey in Figure 7). As already noticed in Remark 4.17, a paired calibration for 
the partition (A 2 , A 2 , A 3 ) is given by the constant vectorfields <pi := ^(-^,3), 
(f>2 := |(— 2v / 3,0), and 4> 3 := — 3), but we can take solutions tpi of (4.13) 

such that < 2e on Q, independently of the value of r. More precisely, we 
take the solution ipi of (4.13) which takes the value on the transversal set Mi 
described in the Figure 7 for i = 1 (the construction is similar for i = 2). 



CI (£-neighborhood of Su) 
integral curves of <j> 2 -<p\ 

transversal set M l 

singular set Su 

Figure 7 

Therefore, if a i+1 — > 2\f2e, the construction of Paragraph 4.16 yields a cali- 
bration of u on Q, x R. Moreover this calibration can be extended to a calibration 
on an e-neighbourhood of the complete graph of u over the entire B(0, r) by set- 
ting it equal to where it is not yet defined. Clearly, this requires the slabs U{ 
are contained in B(0, r) x (a^ + e, aj + i — e) for i = 1, 2, that is 




H+l 



>2e + 2V2e fori = 1,2. (4.18) 



This shows that u is a Dirichlet JZ-minimizer of F when U is an £-neighbourhood 
of the complete graph of u within B(0, r) x R and e satisfies (4.18). As expected, 
e does not depend on the size of the domain, but only on the relative distances 
of the values a,-. 



5. Applications to the complete Mumford-Shah 
functional 



In this section we focus on minimizers of the complete Mumford-Shah func- 
tional F defined in (1.1) for a, (3 > 0, and calibrations will always be intended as 
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in the sense of Theorem 3.4 and Definition 3.5. All the following examples are in 
dimension n. 

5.1. Solutions of the Neumann problem. - If we restrict F to functions 
of class W 1 ' 2 , we obtain the strictly convex and coercive functional Jq[|Vu| 2 + 
(3(u — g) 2 } dx, and its unique minimizer u is the solution of the Neumann problem 

Au = (3{u — q) on fL 

V ' (5.1) 
d v u = on d£l, 

where d v denotes the normal derivative. Thus it is natural to ask under which 
assumptions (on g and (3) u is also a minimizer of F on SBV(Q). This question 
is akin to the minimality of harmonic functions for F discussed in Paragraph 
4.10, and following the same ideas we can construct an absolute calibration for u 
provided that u satisfy condition (5.3) below (cf. also Paragraph 5.4). 

More precisely, we assume that u is of class C 1 up to the boundary (this is 
always satisfied if dfl is of class C 1,e for some e > 0), we denote the infimum and 
the supremum of g by m and M respectively, and set 

A:={(x,t)enxll: ^<f<^}, 



and (cf. (4.9)) 

r(l , t):= ( 2V »W * (5.2) 

1 otherwise. 



Note that, by the maximum principle, m < u < M on £1, so that Tu is contained 
in A. Independently of the choice of <j) t , we can already see that assumption (b) 
of Lemma 3.2 is satisfied if (cf. (4.8)) 

osc g ■ sup | Vu\ < a , (5.3) 

while assumption (b') is trivially satisfied, and 4> has vanishing normal component 
on dfl x R (because of (5.1)). Thus it remains to choose cf) 1 so that (a) and (a') 
hold, and 4> is approximately regular and divergence-free. 

Assumption (a') sets 0* equal to |Vu| 2 — (3(u — g) 2 on the graph of u., while 
requiring that (f> is divergence-free in the interior of A yields 

dttf = -div x <f) x = -2Au = -2/3(u - g) . 

Integrating in t we obtain that is given in A by 

<? = |Vu| 2 - [3(u - g) 2 - 2(3(u - g)(t - u) 
= \Vu\ 2 -f3(t-g) 2 + P(t-u) 2 . 
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Therefore assumption (a) of Lemma 3.2, namely 4> l > \Vu\ 2 — (3(t — g) 2 , is clearly 
satisfied in A. Moreover cf> is approximately continuous in A (this is trivial for the 
vectorfield (4> x ,0), which is continuous, and follows from Remark 2.3 for (0, 0*)). 

Moreover, cf> is divergence- free in the complement of A if we impose that d t $ = 
0, that is, 4> l depends only on x, while the compatibility condition (2.5), which is 
required in order to have that <j) is divergence-free on the entire Q x R (cf. Remark 
2.8(b)), yields 

r^(H±M_ 5 ) 2 + ^(^) 2 fort>M±M, 
" \-/3(^ - g) 2 + (3(^) 2 forf<^. 

Finally one can easily check that 4> is approximately regular also outside A (see 
Remark 2.3) and satisfies condition (a) of Lemma 3.2 as well. Therefore 4> is an 
absolute calibration for u, provided that u satisfies (5.3). 

Remark 5.2. - If inequality (5.3) is strict, the calibration constructed in the 
previous paragraph satisfies the strict inequality in assumption (b) of Lemma 3.2, 
and by Remark 3.6 this proves that the solution u of (5.1) is the unique absolute 
minimizer of F. 

Remark 5.3. - The weaker inequality 

oscu • sup |Vu| < a , (5.4) 

which follows from (5.1) and (5.3) by the maximum principle, is not enough to 
guarantee the minimality of a solution u of (5.1), not even the Dirichlet minimality. 
Take indeed n := 1, Q. := (—a, a), and g(x) := 1 for x > g(x) := —1 for x < 0. 
Then the solution to (5.1) can be computed explicitly: 

where 7 := \f]5a. Now we fix a so that (5.4) is satisfied, and precisely 

1 



a := oscu • sup \ Vu\ = 2\f]3 tanh 



7 



1 - 



cosh 7- 



and take the comparison function v (x) := [l — 1/ cosh 7] g(x); v has the same 
boundary values as u, and a tedious but straightforward computation gives 



F(v) = 2a//3 tanh 7 - 2 ^ (simi7 - 7) < 2v7? tanh 7 = F( 

cosh 7 



Therefore u satisfies condition (5.4) but is not a Dirichlet minimizer of F. 

5.4. Solution of the Neumann problem for large (3. - The construc- 
tion in Paragraph 5.1 shows that the solution of the Neumann problem (5.1) is 
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an absolute minimizer of F provided that (5.3) holds. However, this condition 
is far from being necessary. In particular, for large values of the penalization 
parameter (3, the absolute minimizer u of F is close to g, and therefore we expect 
that discontinuities should not be energetically convenient, at least for sufficiently 
regular g, and the solution u of (5.1) should be the unique absolute minimizer 
of F. 

We prove this fact by calibration under the assumption that £1 has boundary 
of class C 2 , g is of class W 2 ' p for some p > n, and (3 is larger than a certain (3 , 
specified in (5.15). Under these assumptions, g belongs to C x,1 {0) for 7 := l—n/p, 
and u belongs to C 3 ' 7 (ft) n C M (H) n W 2 ' P (Q) for every 5 £ (0, 1) by the standard 
regularity theory for Neumann problems (see, e.g., [39], Theorems 3.5, 3.16, and 
3.17). 

Fix a positive constant 5 (to be properly chosen later), and take a smooth 
function a : R — > [0,1], with support included in [—2(5,25] and identically equal 
to 1 in [—5,5], so that |<r| < 2/5 (and then ||cr||i < 45 and ||<t||oo < 2/5). Set 

4> x (x, t) := 2a(t - u(x)) Vu{x) . (5.5) 

To simplify the notation, in the following we simply write a and Vu instead of 
a(t — u(x)) and Vu(x) (this must be kept into account when deriving), so that 
(5.5) becomes simply <fi x = 2a Vu. 

It follows from (5.1) and (5.5) that <j) has vanishing normal component at the 
boundary of Q x R, and 4> x = 2Vu on the graph of u. Assumption (a') in Lemma 
3.2 prescribes the value of on the graph of u, and precisely 

<j>*{x,u) := |Vu| 2 -(3{u-g) 2 for all x G H. (5.6) 

We impose now that 4> is divergence-free, which reduces to 

dttf = -dvv x (f> x = -2aAu + 2a\Vu\ 2 

= -2f3a{u- g) + 2a\Vu\ 2 . (5.7) 

Identities (5.7) and (5.6) together determine 4> l everywhere. 

Now we want to verify that assumption (a) of Lemma 3.2 holds, that is, 

^>\m 2 -m-9? ■ (5.8) 

Since the equality holds by construction on the graph of u, the full inequality is 
implied by the following inequalities on the derivatives of both sides of (5.8) with 
respect to t: 

dtft >\cp x d t (j) x -2l3{t- g) fori>u, 
<9 t 0* < \(j) x d t (t) x -2(3{t-g) iort<u. ^'^ 
Let us consider the first inequality: by (5.5) and (5.7) it becomes 

-2/3ct (u-g) + 2d \Vu\ 2 > 2aa | Vu| 2 - 2(3{t - g) , 
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which is equivalent to 

f3[(t -g)- a{u - g)\ > &(a - l)\Vu\ 2 . (5.10) 

When jj < K 11 + 15 we have a = 1, and then (5.10) becomes t — u > 0, which is 
obviously true. When t > u + 5, we have (t — g) — a(u — g) > 5 — \\u — g\\oo an d 
|cr(cr — 1)| < 2/5, and then (5.10) is implied by (3(5- \\u- gW^) > f||Vu||^. This 
inequality can be rewritten as 

<5 2 -^h-5||oo-|||Vn||^ >0 , 

and is satisfied for 

[2 

S > \\u- plloo + y - || Vti||oo ■ (5.11) 

One checks in the same way that (5.11) implies the second inequality in (5.9) too. 
In other words, assumption (a) of Lemma 3.2 holds if (5.11) holds. 

Assumption (b') of Lemma 3.2 is trivially satisfied because Su is empty, while 
(5.5) and the estimate \\<r\\i < 4(5 imply that assumption (b) of Lemma 3.2 is 
satisfied with strict inequality if S^HVuHoo < a, that is, 

<5< „ "„ . (5.12) 
8||V«||oo 

Finally, we can find 5 that satisfies both (5.11) and (5.12) if 

HVulU^ || u - gWce + v^||V«||oo) < ^VP , (5-13) 

and by Theorem 3.4 and Remark 3.6 we conclude that, if (5.13) is satisfied, then 
u is the unique absolute minimizer of F. 

Thus it remains to show that (5.13) holds for f3 large enough. Note that u, 
being a solution of (5.1), depends on (3, and there exist positive constants K and 
j3 (depending on but not on g and j3) such that for every (3 > j3 there holds 

w^p ■ (5-14) 

This estimate can be derived, for instance, from the interpolation inequality 
(3.1.59) of Theorem 3.1.22 in [28] (one has to replace A, u, A, and B with (3, 
u — g, A, and d u respectively, and recall that Au = f3(u — g)). 
Estimate (5.14) shows that (5.13) holds for 

(3>(3 := max{p, 2 7 a" 2 K 4 || Vg\\^, P } . (5.15) 



5.5. Characteristic functions of regular sets. - If g := l E is the 
characteristic function of a sufficiently regular compact subset E of Cl, then it is 
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natural to conjecture that for large values of (3 the minimizer of F is g itself. We 
prove this statement by calibration under the assumption that the boundary of E 
is of class C 1 ' 1 (cf. also Remark 5.6 below) and (3 > (3q, with (5q defined by (5.22). 
Under these assumptions we also prove the uniqueness of the minimizer. 

As in the previous paragraph, we first construct 4> x . To this end, we take a 
Lipschitz vectorfield tp on Q which agrees on dE with the inner normal of dE, 
is supported on a neighbourhood of dE which is relatively compact in Q, and 
satisfies < 1 everywhere. For instance, we can use that dE is locally a graph, 
which yields a trivial extension of the normal vectorfield on a small neighbourhood 
of each point, and then use a partition of unity to paste together these different 
extensions. Now we set 

(f> x (x,t) := a(t) %j){x) for all x £ U, t £ R, (5.16) 

where a : R — > [0, 2a] is a function of class C , supported in [0, 1], with integral 
equal to a, and such that \&(t)\ < 16a for t £ [0, 1], a(t) := t 2 for t G [0, 1/8], 
alt) ■= (1-t) 2 for te [7/8,1]. 

We see that, independently of the choice of the vectorfield <j) has vanishing 
normal component at the boundary of Q, and satisfies assumptions (b) and (b') 
of Lemma 3.2. Since 4> x vanishes for t = and for t = 1, and therefore on the 
graph of g, requiring that (f> satisfies assumption (a') yields 

<j>*(x,g{x)) :=0 for T-a.e. i e fi, (5.17) 

while requiring that </> is divergence-free yields (cf. (5.16)) 

8^ = -fcv x (j) x = -crdiv^V • (5.18) 

Conditions (5.17) and (5.18) together determine 0*. 

Note that <j) is approximately regular: this is trivial for the vectorfield (4> x ,0), 
which is continuous, and follows from Remark 2.3 for (0,^*) (even though (f) 1 is 
discontinuous on dE x R and where div^?/; is). 

To show that <j) is an absolute calibration it remains thus to verify assumption 
(a) of Lemma 3.2, namely 

v > i\r\ 2 - Pit - g? = - m - 9? ■ (5.i9) 

Since the equality holds by construction on the graph of g (cf. (5.17)), it is enough 
that <9 t </>* satisfies the inequality 

dttf := -o-div^ > \aa\ijj\ 2 - 2(3{t - g) for t>g(x), (5.20) 

and the opposite inequality for t < g(x). Inequality (5.20) is clearly satisfied for 
t > 1, since a(t) = 0. If g(x) = and < t < 1, (5.20) is implied by 

~cr||div x V||oo > ^\&\ - 2(3t . (5.21) 



36 



G. Alberti, G. Bouchitte, G. Dal Maso 



In turn, (5.21) reduces for < t < 1/8 to 

-^Hdiv^Hco > t 3 -20t , 

which is satisfied for (3 > Hdiv^lloo + j^g, while, for 1/8 < t < 1, (5.21) follows 
from 



which is satisfied for (3 > SaHdiv^Vlloo + 64a 2 . Therefore (5.20) holds for 



The same condition implies also the opposite inequality for t < g{x). This con- 
cludes the proof that cj) calibrates g. 

To prove that g is the unique minimizer of F, we first notice that the strict 
inequality in (5.20) implies the strict inequality in (5.19) for t > g(x), and of course 
we have the strict inequality for t < g(x), too. In other words, the inequality in 
assumption (a) of Lemma 3.2 is strict for all t ^ g{x). Now, if u is another 
minimizer, <fi must calibrate u, too (cf. Remark 3.6), and in particular it must 
satisfy assumption (a') of Lemma 3.2 for u, which means that the inequality in 
assumption (a) is an equality for t = u(x). Therefore we conclude that u(x) = g(x) 
for £ n -a.e. x in Q. 

Remark 5.6. - If g := 1e is the characteristic function of a set E relatively 
compact in VL and u := g minimizes F, then the set E minimizes in particular 
F(A) := F(1 A ) = aH n_1 (<9*,4) + 0\AAE\ among all sets A with finite perimeter 
in f2. Hence the regularity theory for minimal perimeters yields that that, in 
dimension n < 7, E must be of class C 1 ' 7 for every 7 < 1, while in dimension 
two it must be of class C 1 ' 1 (see, e.g., [5], Theorem 4.7.4). Thus the regularity 
on g required in the previous paragraph is optimal in dimension two, and close to 
optimal for 3 < n < 7. 

Remark 5.7. - If g is the characteristic function of a set E of class C 1 ' 1 
which is not relatively compact in Q, then the result of Paragraph 5.5 can be 
generalized as follows: u := g is a minimizer of F for large values of (3 provided 
that BE is orthogonal to dQ, (which is assumed to be sufficiently smooth). An 
absolute calibration can be constructed as in the previous paragraph, one has only 
to choose tp so that it is tangent to the boundary of Q. 

Notice that this orthogonality requirement is necessary: indeed it is easily 
proved that, given a minimizer of the functional aH n ~ 1 (Q n d*A) + f3\AAE\ 
among all finite perimeter sets A in Q, its boundary is orthogonal to <9H 

We conclude this section with some remarks on the gradient flow associated 
with the (homogeneous) Mumford-Shah functional. 

5.8. Gradient flow for the Mumford-Shah functional. - A gradient 
flow for F with respect to the L 2 -metric can be defined in a variational way by 



-2a||div a! V||oo > 16a 2 - \f3 




(5.22) 



The calibration method for the Mumford-Shah functional 



37 



time discretization, following the minimizing movements approach developed in 
[19], [4], [27]. Given an initial datum uq 6 L 2 (Q.) and a discretization step 5 > 0, 
we set m<5 5 o : = and define inductively u§j for j = 1, 2, ... as any minimizer of 

F (u) + ^ Jju-usj-ifdx (5.23) 

among all functions u in SBV(Q) - with or without prescribed boundary values, 
according to the boundary condition (Dirichlet or Neumann) imposed on the 
flow. Then we define u$ : Q x [0, +oo) — > R by u$(x, t) := usj for t := jS, and by 
linear interpolation for t G (j<5, (j + 1)<5), and call gradient flow with initial datum 
u(x, 0) = uo(x) any possible limit of us as 6 — > along any sequence. Note that 
the flow may be not unique, as even u§ is not uniquely defined. 

Remark 5.9. - If uq belongs to W 1,2 (£l), and the minimization of (5.23) is 
restricted a priori to the functions u in W 1,2 (Q), then F (u) is just the usual 
Dirichlet integral, and it can be proved (cf. [4], Example 2.1) that the previ- 
ous construction yields a unique flow which agrees with the solution of the heat 
equation 

d t u = Au on Q x (0, +oo) 

with initial datum u(x, 0) = uq(x) and boundary conditions - Neumann or Dirich- 
let - according to the boundary conditions imposed in the minimization of (5.23). 

The previous remark and the result of Paragraph 5.4 suggest that for a smooth 
initial datum no, the gradient flow associated with Fq is just the solution of the 
heat equation. To prove this, however, we need some additional information on 
the minima of F for large (3. 

5.10. Improved estimates on the solution of (5.1). - Under the regu- 
larity assumptions on £1 and g given in Paragraph 5.4, if Ag 6 L°°(Q) and d v g = 
on dfl, then the solution u to the Neumann problem (5.1) satisfies 

||Au||oo < HAsHoo , (5.24) 

and an improved version of estimate (5.14): 

/% " Slloo + UVulloo < tfllA^loo , (5.25) 

with K depending on Q,, but not on g and (5. In particular condition (5.13) of 
Paragraph 5.4 holds for 

P > A> := max{l, 2 7 a~ 2 K 4 \\ A ff ||^} , (5.26) 

and in that case u is the unique absolute minimizer of F. 
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To prove (5.24) and (5.25), we first notice that the function v := g + e is a 
super-solution of (5.1) as long as e > /3~ 1 \\Ag\\ 00 , in the sense that 

At; < (3{v - g) on fi, 
d v v < on <9H 

Thus u is (a.e.) smaller than g + e on Q. Similarly, g — e is a sub-solution, and 
then 

Ih-Slloo < ^||A 5 |U , (5.27) 

which, in view of (5.1), implies (5.24). 

Now, u solves the equation Am = / with Neumann boundary conditions and 
/ := (5{u — g), and well-known estimates (cf. [39], Theorem 3.16) give ||Vu||oo < 
^"ll/lloo f° r a suitable constant K depending on Q, but not on /. Together with 
(5.27), this implies (5.25). 

5.11. Solution of the Dirichlet problem for large (3. - Assume 
that SI and g satisfy the regularity assumptions of Paragraph 5.4, Ag belongs 
to L°°(Q), and go is a function in W 2,P (Q), and consider the solution u to the 
Dirichlet problem 

Au = 8{u — q) on Q, 

V ' (5.28) 
= go on d£l. 

Then u belongs to C 3 'T(fi) n C 1 ' 5 ^) n W 2 ' P (Q) for every 5 £ (0, 1) (see, e.g., [39], 
Theorems 3.5, 3.16, and 3.17). We claim that if go = g on d£l and (5 is sufficiently 
large, then u is the unique Dirichlet minimizer of F with boundary value go. 

This claim can be proved by the same calibration constructed in Paragraph 
5.4, provided that estimate (5.14) is suitably replaced. To this end, we notice 
that v := g + e is a super-solution of (5.28) for e > /3 _1 ||Agi|| 00 , in the sense that 

Av < @(v - g) on f2, 
v > go on dft 

(we use here that g = go on dfl), and, similarly, g — e is a sub-solution. As in the 
previous paragraph, we deduce that (3\\u — g\\ < WAgH^, and hence (cf. (5.24)) 

||Au||oo < ||A 5 || 

Let us consider now the function w := u — go', since it solves Aw = f — Ago 
with Dirichlet boundary conditions w = 0on dfl, and / := (5{u — g), well-known 
estimates for solutions of Dirichlet problems (see, e.g., [39], Theorem 3.16) give 
|| Vty||oo < i^dl/Hoo + ||V<7o||c >^)- Together with the estimate on \\u — g\\oo, this 
implies (cf. (5.25)) 

0\\u -g\\oo + II V«||oo < ^(11 Aslloo + ||V ff0 || c o, 7 ) , 
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with a possibly different K. In particular, condition (5.13) of Paragraph 5.4 is 
satisfied for 

(3 > fa := max{l, 2 7 a - 2 K 4 (\\ A^U + ||V<7 ||co-) 4 } , (5.29) 

and in that case u is the unique Dirichlet minimizer of F with boundary values 
u = g - 

5.12. Gradient flow with smooth initial datum. - Assume that has 
boundary of class C 2 , Uq G W 2 ' p (Q) for some p > n, Auq G L°°(Q), and d v uo = 
on dVt. Then the gradient flow for F with Neumann boundary conditions and 
initial datum u(x,0) = uq(x) constructed in Paragraph 5.8 is unique, and agrees 
with the solution of the heat equation. 

In virtue of Remark 5.9, this claim is an immediate consequence of the follow- 
ing fact: when 



8 <8 



max{l, 2 7 a- 2 K 4 ||An ||^} 



i 



then, for every integer j, every minimizer of (5.23) belongs to W 1,2 (Q). In other 
words, the solution of the Neumann problem (5.1) with (3 := 1/8 and g := usj-i 
is the unique minimizer of (5.23). 

To prove this fact, it suffices to verify that the assumptions of Paragraph 5.10 
are satisfied for every j, and precisely: usj-i G W 2,P (Q), Aus^-i G L°°(Q), 
d u u$j-i = on <9r2, and inequality (5.26) holds with (3 := 1/8 and g := ugj—i- 
The last requirement follows from the choice of 8 and the chain of inequalities 

||A« ||oo : = ||A-U 5j0 ||oo > IIAu^iHoo > || Ati 5)2 ||oo > ••• , 

which are implied by (5.24). The W 2,p regularity of usj follows from the corre- 
sponding regularity of usj-i, as remarked at the beginning of Paragraph 5.4. 

Remark 5.13. - The conclusion of the previous paragraph also holds for the 
gradient flow with Dirichlet boundary conditions. More precisely, if Q has bound- 
ary of class C 2 and uq is of class W 2 ' P (Q), with Auo G L°°(r2), then the gradient 
flow for F with initial datum u(x, 0) = uq(x) on Q and Dirichlet boundary con- 
dition u(x,t) = uq(x) on dU x [0, +oo) constructed in Paragraph 5.8 is unique, 
and agrees with the solution of the heat equation (with same initial datum and 
boundary conditions). The proof is essentially the same as in the Neumann case, 
and relies on the estimates given Paragraph 5.11 (in particular one has to apply 
(5.28) with g = u S j-i and g = u ). 

Remark 5.14. - If we drop the assumption 3 u uq = on 90 in Paragraph 
5.12, the proof breaks down because we can no longer estimate || Augj-i by 
||A«o||oo 5 but we do not know if the conclusion on the gradient flow still holds. 
A similar problem occurs with the statement in Remark 5.13 if we replace the 
Dirichlet boundary condition u(x,t) = Uq{x) on dQ with a different one. 
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Remark 5.15. - If Uq := 1e is the characteristic function of a compact subset 
of Q. with boundary of class C 1 ' 1 , then the gradient flow for F with Neumann 
(or Dirichlet) boundary conditions constructed in Paragraph 5.8 is unique, and 
is just given by u(x,t) := Uo(x) on ft x [0, +oo). This follows immediately from 
Paragraph 5.5. 

Remark 5.16. — The conclusions of Paragraph 5.12 and Remark 5.15 sup- 
port the following general conjecture: if the initial datum uq is smooth out of a 
smooth singular set Su, which is compact in Q, then, at least for small times, 
the associated gradient flow should be unique and should keep the singular set of 
u(-,t) equal to Suq, while the function u evolves in Q \ Su according to the heat 
equation (with Neumann conditions on Su). This conjecture has been proved 
in the one-dimensional case in [26] (with a slightly different definition of gradient 
flow for F ). The general case will be studied in [35] using the calibration method. 

6. Appendix 

In this section we prove some technical lemmas stated in Section 2. We follow 
the notation of Section 2. 

Lemma 6.1. - Let Q be an open subset of R n whose boundary is the graph 
of a Lipschitz function f : R n_1 — > R, and let 4> be a bounded vectorfield on 
which has bounded support and satisfies condition (2.4) with M := dQ. Then 
there exists a sequence of vectors yj such that yj — > 0, dfl + yj C for every j, 
and 



PROOF. Let S be the set of all vectors y e R n such that dfl + y C ft. For 
every r > 0, let S r := S PI B(0,r), and consider the double integral 



If we invert the order of integration, condition (2.4) means that the inner integral 
(over S r ) tends to as r — > for H n ~ 1 -&.e. x G dfl. Then (6.2) converges to by 
the dominated convergence theorem (recall that (f> is bounded and has bounded 
support). 

As d£l is a graph of a Lipschitz function, the set S contains an open cone with 
vertex in 0. Then the measure of S r is larger than ar n for some fixed a > 0, and 
therefore we can choose y r 6 S r so that the value of the inner integral in (6.2) is 
smaller than the double integral divided by a, and then converges to as r — ► 0. 

In other words, <f>(x + y r ) ■ u dn (x) converge to 4>{x) ■ v QQ (x) in L 1 (9f^, H n_1 ), 
and then it suffices to choose a subsequence yj which yields pointwise convergence 
for H^-a.e. x £ dfl. □ 



lim 4>(x + yj) ■ u dn (x) = 4>{x) ■ v dn (x) forTL n -a.e. x 6 dQ. 



(6.1) 




(6.2) 
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PROOF of Lemma 2.4. - (Sketch) We divide the proof in several steps. 

Step 1. Assume that belong to C,!(R n ,R n ). In this case formula (2.6) is 
well-known - see, e.g., [25], Theorem 2.10, or [6], formula (3.87) in Theorem 3.87. 

Step 2. Assume that 4> is an approximately regular vectorfield on R n with 
compact support and that divcj) 6 L°°(R n ). Let ip s (x) := e~ n ip(x/e) be a stan- 
dard radially symmetric regularizing kernel of class C£°, and take <p £ := (f) * ip e . 
Thus formula (2.6) holds for each <fi £ by Step 1, and it only remains to check that 
we can pass to the limit as e — > 0. The convergence of the first integral in the 
right-hand side of (2.6) follows from the fact that the functions div<j) E are bounded 
in L°° and converge to div^ a.e. in Cl. Since </> is approximately regular, the maps 
4> e ■ v M converge to <fi ■ v M H n ~ 1 -a.e. on any Lipschitz surface M, and then also 
on any rectifiable set M. In particular this implies the convergence of the second 
integral in the right-hand side of (2.6). The same argument also applies to the 
left-hand side, provided that we use the coarea formula (cf. [6], Theorem 3.40, or 
[23], Theorem 4.5.9(13)) to re-write that integral as 



where M t is the measure theoretic boundary in Cl of the sublevel {u < t}. 

Step 3. If ^ is a compactly supported, approximately regular vectorfield on a 
neighbourhood of Cl with divcp in L°°, we reduce to Step 2 using a suitable cut-off 
function. 

Step 4. Assume that Cl is the subgraph of a Lipschitz function / : R n_1 — > R, 
and is a compactly supported, approximately regular vectorfield on Cl with 
div(f) € L°°(Cl). We take a sequence of vectors yj as in Lemma 6.1, and set 
4>j(x) := <p(x + yj), Uj(x) := u(x + yj). By Step 3, formula (2.6) holds with <f> and 
u replaced by <pj and Uj , and it remains to check that we can pass to the limit as 
j — > +oo. The convergence is immediate for all integrals in (2.6) but the last one. 
In this case, it suffices to notice that the functions 4>j ' v dn are uniformly bounded 
and converge to <fi ■ v gn H n_1 -a.e. on dQ (by the choice of the vectors yj), while 
the traces of Uj on dfl converge to the trace of u in L l (dVL,H n ~ l ) (because the 
functions Uj converge to u in variation, or, alternatively, because the L 1 -norm of 
the difference of the traces is controlled, up to a constant which does not depend 
on j, by |Z>u|(f2 \ (CI — yj)), which clearly tends to zero). 

Step 5. To prove the general case, we use a locally finite partition of unity 
consisting of compactly supported smooth functions to reduce to Step 4. □ 

PROOF of Lemma 2.6. - We first prove that div^ = / on Ci \ So- Since the 
problem is local, it is enough to show that div0 = / on every bounded open set 
U with U C Cl \ So and such that U\S\ has two connected components U + and 
U~ with Lipschitz boundary. As (p is approximately regular on U^, we can apply 
formula (2.6) with Cl replaced by and u 6 C^°(U): as the integrals on U n Si 
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cancel out, we are left with J v 4> • Vu dx = — f v fu dx. By the arbitrariness of u, 
we deduce that div0 = / on U. 

We prove now that div0 = / on Cl. Since H n ~ 1 (So) = 0, the (1, 1) -capacity 
of So is zero (see [22], Section 5.6.3), and therefore there exists a sequence of 
functions Gj in C°°(Q) such that < Uj < 1 in Q, and Uj = in a neighbourhood 
of Sq, and aj — > 1 strongly in W 1 ' 1 (f2). 

Take now an arbitrary function u G C^°(fl). Then the functions UjU be- 
long to C^°(Q \ So) and, since divc/> = / on O \ So, we have J n • V(aju) dx = 
— / ■ (<7j"it) dx. Moreover the functions aju converge to u strongly in W ' (ft), 
and therefore f n 4> ■ Vudx = — / n fudx, which concludes the proof. □ 

PROOF of Lemma 2.10. - By a monotone class argument it is enough to 
prove (2.9) for <j) of the form <p(x,t) := p(t)ip(x), with p : R — > R and -0 = 
(ip x ,ip t ) : f2 — » R n x R of class C£°. Let <r be the primitive of p vanishing at — oo. 
Then we have 



4> ■ Dl u = — div0 l u dx 

fixR JttxR, 



(pdiv x ip x + pV*) (it cfcc 
= - f [a(u) div x ip x + p{u)i) 1 ] dx . (6.3) 

As u belongs to SBV(Q), the chain-rule for BV- functions (see, e.g., [6], Theorem 
3.96) gives 

D[a(u)} = p(u)Vu-C n + [a(u + )-a(u-)]iy u -n n - 1 ^Su . 
Therefore (6.3) implies 

f cp ■ Dl u = [ [p{u)^ x ■ Vu - p(u)iff*] dx 

+ [ [a{u + )-a{u-)}i> x -u u dH n - 1 , 

JSu 



which, together with (2.8), gives (2.9) in the case <fi(x,t) := p(t) tp(x). □ 

PROOF of Lemma 2.12. - We set w := l u — l v on x R. Then w belongs 
to BV (fi x R) and Dw = v Vu -H n LLu - v Tv -H n L r«. 

Let us consider the inner trace of u> on First of all we decompose dU as 
the disjoint union of (Q x R) n dU and (dfl x R) n dU. For every C°° vectorfield 
ip on x R with compact support we apply formula (2.6) of Lemma 2.4 with O 
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and 4> replaced by U and ip, respectively, and we obtain 

— w div-0 dx = I ip ■ Dw + w ip ■ u QU dH n 

Ju JU JdU 

= [ i>-v Vu dH n - \ V • y Tv dU n (6.4) 

JFunU JTvCtfJ 



+ J wtp ■ v QV dTi r 
hu 



On the other hand, by the definition of distributional derivative we have also 



/ w div^ dx = ip ■ Dw 

JnxR, isixR 



f V • v Tu dH n - [ ^ • v Vv dU n . 

JFu JFv 



(6.5) 



Due to the particular structure of U and the assumption on the complete graphs 
of u and v, the function w vanishes a.e. on (fi x R) \ U. This fact, together with 
(6.4) and (6.5), implies that the inner trace of w on (17 x R) n dll satisfies 



w v 



dU 



l Tu v Vu - \ Vv v Vv TT-a.e. on (SI x R) n dU 



Therefore w belongs to ^((Q x R) n dU,H n ) and 

f w^-v m dU n = [ <j> ■ v Vu dU n - [ </>■ v Tv dH n . (6.6) 

J(fix~R.)ndu JFundU JFvDdU 

Now, the trace of w on dSl x R is the difference of the characteristic functions 
of the traces of u and v on dSl, and therefore it belongs to L 1 (90 x R, H n ) and 
vanishes H n -a.e. on (d£l x R) \ dU. As uqu = (^an,0) on (dCl x R) n dU, this 
implies 



/ w<j>- v au dTT = w<jf ■ v dn dTL 

J(dnxR,)ndU ' JdnxR '" (6 7) 

4> x (x,t)dt 



v^dH"- 1 . 

Since the inner trace of w on dU belongs to L 1 (dU,H n ), we apply formula (2.6) 
of Lemma 2.4 with Q and u replaced by U and w, respectively, and get 

/ cf*-v ru dH n - [ <f> ■ v Vv dH n = - f w<\>-v m dH n . (6.8) 

JFunU JFvnU JdU 

Identity (2.10) follows now from (6.6), (6.7), and (6.8). □ 
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